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Project  Summary 


Modern  requirements  for  high  performance  aircraft  require  advanced  control 
ured  flight  vehicles  where  an  increasing  reliance  is  placed  on  computer  control  .?}  s;  .. 
for  flight  performance.  Such  aircraft  will  be  equipped  with  multiple  sensors  and  actu¬ 
ators  and  employ  centralized  control  implemented  by  a  high  speed  digital  comp; 

For  such  systems,  requirements  for  reliable  multiloop  control  laws  which  are  robust  io 
a  variety  of  changes  in  flight  dynamics  are  prerequisite  for  successful  flight  operation 
Typical  methods  for  the  synthesis  of  high  performance  feedback  control  is  based  m 
nominal  dynamic  model  for  the  aircraft  flight  dynamics  Such  models  are  subject  \  • 
various  limitations  including  the  effect  of  changes  in  flight  conditions  affecting  trim  <  r 
equilibrium  and  consequent  changes  in  the  linear  perturbation  dynamics.  Typical  con¬ 
trol  design  methods  are  based  on  the  linear  dynamics  which  will  be  subject  to  vaii  /c 
-  mrces  of  model  uncertainty  in  flight  control  applications. 

In  this  project  we  have  initiated  studies  of  the  application  of  new  and  advanced 
nethods  for  control  law'  synthesis  for  robust  stabilization  with  respect  to  a  combination 
->f  both  unstructured  model  uncertainty  (arising  from  neglected  or  parasitic  dynamics) 
uul  structured  model  uncertainty  (arising  from  parametric  variations  which  occur  as 
i'ght  conditions  change.)  Our  efforts  have  focused  on  the  characterization  of  a  class  of 
-mnlinear  models  for  longitudinal  dynamics  of  aircraft  in  level  flight  subject  to  changes  in 
tatic  stability.  Such  “relaxed  stability”  aircraft  configurations  are  currently  at  issue  in  a 
vide  variety  of  advanced  designs  including  commercial  transport  and  high  performance 
.arc  raft. 

Our  approach  for  control  design  is  to  employ  H°°  synthesis  methods  for  optimal 
>  .bust  stabilization  for  the  unstructured  model  uncertainty  using  a  computationally 
actable  approach  of  (Hover  and  McFarlane.  The  requirement  for  worst  case  [ i  p  to 
bust,)  design  to  parametric  uncertainty  is  included  using  a  minimax  optimization  criteria, 
oir  studies  have  focused  on  a  class  of  flight  control  models  with  physically  based  para 
•  H-tric  uncertainty.  For  these  models  we  have  the  solution  of  the  minimax  or  worst  case 
ilesign  by  a  straightforward  procedure  which  can  be  readily  combined  with  the  require¬ 
ments  for  robustness  to  parasitic  dynamics  using  the  closed  form  solution  f  the  .pinna: 
robust  stabilization  method  developed  by  Glover  and  McFarlane. 

Robust  design  for  stability  can  be  combined  with  requirements  for  robust  perfor¬ 
mance  using  frequency  dependent  constraints  on  the  system  gains.  In  this  study  we 
highlight  the  extension  of  the  optimization-based  methods  for  robust  design  with  re 
;>ect  to  performance  and  stability.  Examples  are  included  to  illustrate  rhe  methods 


Significantly,  these  methods  allow  frequency  dependent  weighting  (i.e.,  “loop  shaping”) 
to  be  performed  without  specific  regard  to  phase  and  can  be  applied  with  equal  ease  to 
both  single  and  multi-loop  design  problems. 

In  our  modeling  studies  we  have  also  focused  on  the  relation  between  the  basic  state¬ 
ment  of  the  control  design  problem  and  the  nonlinear  aircraft  model  subject  to  physically 
based  parametric  uncertainty.  We  have  identified  a  relationship  between  the  choice  of 
actuation  and  sensing  (i.e.,  inputs  and  outputs),  the  nonlinear  aircraft  dynamics,  and  * 

the  linear  perturbation  models  used  for  control  synthesis  which  highlights  the  limitations 
of  robust  control  synthesis.  In  particular,  we  identify  a  relationship  between  conditions 
for  the  existence  of  static  bifurcations  of  the  equilibria  and  the  location  of  transmission 
zeros  for  the  linear  perturbation  model.  Among  many  questions  which  this  analysis 
can  help  to  address  we  identify — under  realistic  scalings  of  a  generic  relaxed  stability 
aircraft — the  limits  of  linear  control  synthesis  for  parametric  uncertainty. 

We  also  demonstrate  the  computational  simplicity  of  the  proposed  method  for  com¬ 
bined  structured/unstructured  uncertainty  synthesis  for  the  relaxed  static  stability  air¬ 
craft  model.  Comparisons  are  given  with  several  recently  developed  methods  for  robust 
control  design  for  linear  system  models  subject  to  various  combinations  of  parasitic  and 
parametric  model  uncertainty. 
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1.  Introduction  and  Background 


Future  aerospace  vehicles  will  be  expected  to  perform  under  circumstances  in  which  the  dynamics 
of  these  vehicles  are  not  expected  to  be  well  known.  Moreover,  because  they  will  be  required  to 
function  close  to,  or  even  beyond,  open  loop  stability  limits,  it  is  essential  that  their  flight  control 
systems  will  be  able  to  cope  with  modeling  uncertainty. 

The  important  considerations  are  typically  of  two  types,  variation  of  modeled  system  parameters 
and  unmodeled  parasitic  dynamics.  Some  examples  of  the  former  may  be  found  in  the  pitch  axis 
dynamics  of  the  AFTI/F-16  [52],  the  longitudinal  dynamics  of  an  advanced  subsonic  transport 
airplane  [1],  and  the  vertical  dynamics  of  a  helicopter  [2].  Examples  of  the  latter  include  elastic 
structural  deformations  of  the  airframe  [3]  as  well  as  sensor  and  actuator  dynamics  [4].  Parametric 
variations  are  structured  uncertainties ,  whereas  unmodeled  dynamics  are  unstructured  and  typically 
associated  with  time  scale  separation  or  high  frequency  effects. 

Of  special  concern  are  those  situations  in  which  the  aircraft  is  intentionally  designed  to  be  open 
loop  unstable  or  to  have  markedly  reduced  stability  margins.  In  such  cases  parameter  variations 
dramatically  affect  the  dynamics  of  the  vehicle  and  its  flyability.  It  is  necessary  to  equip  such 
aircraft  with  stability  augmenting  feedback  controllers  which  shape  the  handling  qualities  and 
reduce  the  sensitivity  to  parameter  variations.  McRuer  et  al  [4]  provide  a  thoughtful  assessment  of 
the  current  status  and  deficiencies  in  such  systems.  They  note  the  following  points: 

(i)  Controllers  tend  to  be  of  wide  bandwidth,  resulting  in  sensitivity  to  vehicle  structural 
flexure,  actuator  dynamics  and  other  high  frequency  effects;  viz.  unstructured  model 
uncertainty. 

(ii)  While  the  feedback  controller  may  be  designed  to  produce  desired  responses  to  pilot 
commands,  responses  to  external  (atmospheric)  disturbances  may  be  unusual  and 
deleterious. 

(iii)  The  effective  handling  qualities  introduced  by  augmented  feedback  control  systems 
have  not  been  thoroughly  investigated  and  opportunities  for  “task  tailored  dynamics” 
should  be  exploited. 

(iv)  Flying  qualities  may  deteriorate  substantially  near  the  limits 
of  control  effectiveness. 
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The  fact  that  such  compensators  for  control  configured  aircraft  are  of  high  bandwidth  is  repeatedly 
noted  in  the  flight  control  literature  [3].  It  is  also  observed  that  many  approaches  to  design  of 
stabilizing  controllers,  robust  with  respect  to  structured  uncertainty,  often  lead  to  high  gain 
controllers.  However,  Schmitendorf  [5]  illustrates  by  example  that  high  gain  is  not  a  necessary 
element  of  robust  control.  It  is  clear  that  methods  are  needed  which  simultaneously  address  design 
for  both  structured  and  unstructured  model  uncertainty. 

Certain  nonlinear  issues  are  central  to  the  design  of  controllers  for  relaxed  static  stability  aircraft. 
As  already  noted,  augmented  aircraft  performance  may  radically  diminish  near  limits  of  control 
authority.  Methods  of  control  system  design  are  needed  which  explicity  account  for  control 
saturation.  Relaxed  stability  aircraft  operate  near  equilibria  with  small  domains  of  attraction  and, 
consequently,  the  ability  to  recover  from  disturbances  may  be  dramatically  improved  with 
nonlinear  feedback  laws.  An  example  is  given  by  Garrard  and  Jordan  |6j  in  which  the  ability  to 
recover  from  stall  is  significantly  improved  for  a  fighter  aircraft,  operating  at  a  high  angle  of  attack, 
by  the  use  of  a  nonlinear  controller.  Gain  scheduling  in  terms  of  measurable  parameters  has  been 
found  necessary  [1]  to  achieve  desired  performance  over  the  full  flight  envelope  when  linear 
feedback  is  employed.  Exact  linearization  methods  [48]  provide  means  for  designing  nonlinear 
feedback  laws  which  satisfy  these  requirements.  However,  exact  linearization  is  not  always 
compatible  with  control  authority  constraints  and  very  little  is  known  about  the  suitability  of  these 
methods  under  conditions  of  incipient  instability.  Efficient  methods  are  required  for  the  design  of 
nonlinear  controllers  which  articulate  the  tradeoff  between  performance  and  control  restraints. 

In  recent  years,  numerous  investigators  have  considered  a  variety  of  approaches  to  the 
quantification  of  aircraft  handling  qualitites  in  a  way  which  would  be  useful  in  analytical  control 
system  design  [7],  Such  efforts  take  on  special  significance  for  the  design  of  augmented  flight 
control  systems.  Some  general  discussion  is  given  by  McRuer  et  al  [4|  and  additional  insights  may 
be  found  in  Wilhelm  and  Schafranek  [8].  Nevertheless,  many  basic  questions  need  to  be  answered 
if  the  notion  of  systematically  'designing  in'  desired  flying  qualities  via  feedback  augmentation  is 
to  become  a  reality.  For  example,  common  control  architectures  suggested  for  stabilizing  open 
loop  unstable  aircraft  result  in  nonminimum  phase  zeros  in  the  dynamics  as  seen  by  the  pilot.  This 
issue  has  not  been  previously  addressed  although  the  control  of  nonminimum  phase  plants  is 
notoriously  difficult. 

In  this  report,  we  present  preliminary  results  concerning  the  design  of  linear  controllers  for  flight 
control  applications  which  are  robust  with  respect  to  both  structured  and  unstructured  uncertainties. 
Our  approach  centers  around  two  elemental  themes.  First,  we  are  concerned  with  the  origins  of 


2 


parametric  (structured)  uncertainty  in  aircraft  and  in  the  development  of  an  understanding  of  when 
any  such  parametric  variations  are  critical  and  cannot  be  imbedded  in  the  unstructured  uncertainty 
without  unacceptable  conservatism.  Second,  we  propose  a  method  of  worst  case  design  for 
dealing  with  combined  structured  and  unstructured  model  uncertainty  based  on  a  formulation  of  the 
H°°  optimization  problem  due  to  McFarlane  and  Glover  [33]. 

In  Section  2  we  review  control  design  methods  for  unstructured  model  uncertainty,  thereby 
developing  the  required  foundational  background  for  Section  3.  Section  3  provides  a  review  of 
alternative  approaches  to  extending  these  methods  to  accommodate  structured  mouu  uncertainties 
m  combination  with  unstructured  model  uncertainties  and  also  introduces  our  view  of  addressing 
structured  uncertainties  as  a  minimax  optimal  control  problem.  Section  4  consists  of  a  bench. nark 
design  problem:  robust  control  system  design  for  the  longitudinal  dynamics  of  a  relaxed  sialic 
stability  aircraft.  We  introduce  in  a  rudimentary  way  our  ideas  aoout  the  use  ol  uifuicatiou  oiuujms 
as  a  means  of  articulating  the  fundamental  limits  imposed  by  paiameiric  variation*  in  me  auciah's 
inherently  nonlinear  dynamics.  In  this  example  the  minimax  approach  is  shown  to  be  an  method  of 
control  system  design.  In  Section  5  we  summarize  our  conclusions  and  outline  promising 
directions  for  future  research. 

2.  Robust  Control  Design  for  Unstructured  Model  Uncertainty 

Analysis  of  control  system  stability  plays  a  fundamental  role  in  design  tradeoffs  for  feedback 
compensation.  Classical  design  methods  focus  on  a  frequency  domain  description  of  SISO 
systems  to  articulate  the  tradeoff  between  performance  (e.g.  sensitivity  reduction)  and  stability 
margins.  The  use  of  gain/phase  margins  provide  a  standard  quantifiable  description  of  these 
tradeoffs  directly  in  terms  of  frequency  response  data.  From  experience  with  classical  design 
methods  it  has  become  apparent  that  phase  contributions  and  particularly  phase  errors  due  to 
unmodeled  dynamics  play  a  significant  role  in  the  achievable  performance  and  stability  margins.  In 
particular,  phase  errors  can  ultimately  limit  implementable  gains  and  thus  control  bandwidth  with 
its  intrinsic  relation  to  system  performance.  For  MIMO  system  design  problems  there  is  no  simple 
stability  margin  concept  which  can  separately  address  the  change  in  stability  due  to  phase  errors  in 
the  model. 

Recently,  practical  design  for  MIMO  systems  has  focused  on  the  natural  extensions  of  gain  margin 
concepts  to  the  multiloop  case  using  the  spectral  matrix  norm.  Model  uncertainty  for  a  nominal 
plant  model  G(s)  may  be  described  in  an  unstructured  way  by  reference  to  an  absolute  error  (or 
additive  perturbation,  G(s)— >G(s)+Aa(s))  or  a  relative  error  (or  multiplicative  perturbation 
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G(s)— »[I+Am(s)]G(c)).  The  model  error  is  described  by  a  bounding  function  of  an  appropriate 
form  as: 

|M“'S2  S  1  a<“> 

{•  £  I  mN 

Such  measures  of  model  error  are  gross  in  the  sense  that  they  bound  the  model  uncertainty  in  terms 
of  i  conservative  measure  of  the  multiloop  system  gain  as  given  by  the  choice  of  norm. 

For  application  to  design  and  analysis  of  closed  loop  stability  with  uncertain  models  as  above  we 
are  concerned  with  conditions  under  which  the  feedback  will  stabilize  any  one  of  all  possible  plant 
models  satisfying  the  above  bounds.  Consider  the  following  design  paradigm  for  unstructured 
model  uncertainty  115].  Assume  that  a  nominal  plant  model  G(s)  is  used  to  obtain  a  stabilizing 

controller  K(s)  so  that  the  resulting  closed  loop  transfer  function,  H(s)  =  G(s)  (I  +  GK(s))  1 ,  is 

stable.  The  design  is  said  to  be  robust  with  respect  to  stability  if  we  can  show  that  the  fixed  K(s) 
will  also  stabilize  [I  +  Am(s)  ]G(s)  for  any  perturbation  subject  to  the  frequency  dependent  bound. 

Early  developments  in  methods  for  robust  control  design  as  in  [15]  focus  on  requirements  for 
frequency  dependent  shaping  of  the  MIMO  loop  transmission  KG(s)  in  terms  of  the  above  notion 
of  system  gain;  i.e.,  sufficient  conditions  for  robust  stability  are: 

(i)  Vs)  is  a  stable  transfer  function  and 

(ii)  -  "j  f  0  r  0  <  0)  <  oo  jn  region  where  lm»l. 

The  interpretation  of  the  above  conditions  is  intuitively  appealing  in  that  it  generalizes  typical 
engineering  practice  of  loop  shaping  to  multiloop  designs.  In  application  it  may  be  constraining 
for  several  reasons: 

1)  The  characterization  of  system  gain  may  result  in  a  very  conservative  description  of 
model  uncertainty. 

2)  The  restriction  (i)  may  invalidate  results  in  certain  important  applications. 

Initial  attempts  to  generalize  the  notion  of  system  gain  and  thus  permit  analysis  and  design  for  more 
realistic  problems  were  based  on  generalization  of  conic  sector  bounds  for  transfer  functions  [57]. 
In  this  line  of  thinking  the  system  perturbation  was  assumed  to  have  the  form 

A(S)  =  L_1(s)<D(s)R(s) , 
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where  U,R  are  stable  transfer  functions  chosen  to  represent  the  frequency  dependent  model 
uncertainty  in  a  somewhat  more  structured  way  and  <t»  is  a  stable  unknown  transfer  function 
subject  to  the  frequency  response  bound 


IN40)!  -1 
2 

for  all  co.  Application  of  the  small  gain  theorem  [18]  will  justify  various  sufficient  conditions  for 
robust  closed  loop  stability  depending  on  where  we  choose  to  represent  the  model  uncertainty 
within  the  closed  loop  model.  Typical  model  uncertainty  assumptions  and  corresponding 
frequency  dependent  sufficient  conditions  are  summarized  in  Table  1 . 

Table  1:  Conditions  for  robust  stability 

Frequency  dependent  conditions 
_ for  robust  stability _ 

|rGK(I  +  GK)_1L_1](jco  )|  <  1,  V  to 
|[rWG(I  +  M3)  1 L  1  ]  (j  co  )||  <1,  Vto 

3-neglected  plant  modes 

1[rK(»+  GK)  Y1] (j to )ji  <1,  Vto 
G->G+A  2 

In  practical  design  problems,  the  engineer  may  be  forced  to  consider  a  variety  of  mode!  erici 
sources  in  combination,  and  in  many  cases,  certain  forms  of  model  errors  may  be  more  significant 
in  certain  loops.  This  form  of  structured  uncertainty  was  considered  by  Doyle  et  at  [16j.  ilc 
suggested  that  plant  model  uncertainty  could  be  represented  in  a  general  way  by  isolating  the  model 
uncertainty  in  an  additional  external  feedback  loop  as  shown  in  Figure  1.  Thus  the  model 
uncertainty  is  viewed  as  a  parasitic  feedback  from  certain  error  outputs  e  to  certain  disturbance 
inputs  w.  In  this  case  if  we  take  the  transfer  function  model  P,  representing  the  Know  n  dynamic*, 
i.i  the  form, 


Model  Uncertainty 

1  -neglected  sensor  dynamics 
G->  (1+  A)G 

2-neglected  actuator  dynamics 
G->G(I  +  A) 
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the  assumed  form  of  the  uncertain  plant  is  that  of  a  Linear  Fractional  Transformation  (LFT); 

Q[A)  =  Pyu  +  PywA(l  -  PewA)  Veu 

In  this  form  sufficient  conditions  for  robust  stability  can  be  obtained  by  direct  application  of  the 
small  gain  theorem  (using  possible  frequency  dependent  weighting  functions  L,R)  in  terms  of  an 
LFT  of  P  and  K  in  the  form; 

|KPe»v+PeUK(|-Pyu4)‘1pyw)-‘1||  <1 

2 

for  all  real  frequencies. 


Figure  1 :  The  general  structure  of  the  robust  control  problem. 


Then  a  natural  characterization  of  the  significance  of  loop  dependent  modeling  errors  is  to  consider 
the  parasitic  dynamic  model  to  be  structured  in  the  sense  that  the  transfer  function  has  the  form 
A=diag(Ai,...,Am). 

The  association  be  the  spectral  matrix  norm  (understood  as  a  matrix  norm  subordinate  to  the 
Euclidean  vector  norm  '  mite  dimensions)  and  the  maximum  singular  value  makes  contact  with 
the  spectral  properties  of  singular  value  decomposition  [18].  Following  this  line  of  reasoning, 
Doyle  postulated  the  existence  of  a  structured  singular  value  (SSV),  defined  with  respect  to 
particular  block  diagonal  structure  of  the  loop  parasitics.  The  definition  of  the  SSV  proposed  by 

* 
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Doyle  was  chosen  to  emulate  the  spectral  properties  of  singular  values.  In  particular,  the  fact  that 
the  maximum  (resp.  minimum)  singular  values  are  related  to  the  spectral  matrix  norm  as, 

°maiA}  =  N2-  Om|n{A}  =  kT2 

for  A  a  nonsingular  matrix.  Thus,  the  minimum  and  maximum  singular  values  are  useful  in 
bounding  the  spectrum  of  matrices  under  perturbations.  Considerable  effort  has  been  expended  in 
the  search  for  computationally  feasible  algorithms  for  estimating  the  SSV  [58],  however,  as  of  this 
writing  only  fairly  simple  structures  can  be  computed  with  any  accuracy. 

An  important  feature  for  robust  control  design  methods  is  the  incorporation  of  frequency  shaping 
requirements  given  by  any  of  the  various  frequency  dependent  conditions  for  robust  stability 
described  above.  Such  considerations  have  been  implicit  in  design  methods  dating  back  to  early 
!  970's  [28  j.  One  popular  method  provides  an  extension  to  the  state  space  constructions  of  Linear 
Quadratic  Gaussian  (LQG)  optimal  control  problem  by  the  incorporation  of  Loop  Transfer 
Recovery  (LTR)  [19].  We  remark  that  the  limitations  described  above  all  apply  to  LQG/LTR  type 
methods.  The  popularity  of  this  method  derives  from  the  generality  of  approach  to  MIMO  design 
and  the  computational  support  available  from  various  standard  numerical  algorithms  available  for 
its  implementation. 

In  the  Phase  I  study  we  have  investigated  application  of  several  new  results  in  analysis  and  design 
of  robust  control  systems  which  address  the  above  limitations.  First,  we  employ  a  new  type  of 
stability  margin  for  multiloop  feedback  which  quantifies  closed  loop  internal  stability  in  a  perfectly 
general  way  and  embraces  possibly  unstable  model  uncertainty.  The  stability  margin  can  be 
alternately  described  from  a  geometric  viewpoint  or  as  the  supremum  over  frequency  of  the  gain  of 
a  transfer  function  specially  constructed  from  normalized  coprime  factorizations  for  the  plant  and 
the  compensator.  Furthermore,  recent  results  of  Glover  and  McFarlane  provide  a  synthesis 
procedure  for  obtaining  a  controller  which  achieves  (if  possible)  a  prespecified  stability  margin  of 
this  type.  Their  methods  provide  a  simple  computational  scheme  for  determining  the  maximum 
possible  stabilitv  margin  attainable  for  a  given  plant  model.  To  illustrate  the  application  of  these 
results  we  have  focused  on  a  physically  motivated  flight  control  design  example  which  is  discussed 
in  the  next  section. 

2.1.  Geometric  Stability  Margin  and  Normalized  Coprime  Factorization 

The  geometric  view  of  feedback  considers  the  general  MIMO  feedback  equations  with  G(s)  a  pxm 
plant  transfer  function  and  K(s)  mxp; 
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y(s)  =  G(s)u(s) 
u(s)  =  -  K(s)y(s) 


(2.1) 


in  the  form 


l«  Jly(s)J=a  aia) 

In  this  form  we  focus  attention  on  the  relative  orientation  of  two  abstract  objects; 
Gs=ker[G(s),Ip],  Ks=ker[Im,K(s)],  which  can  be  viewed  for  any  complex  s  as  a  pair  of 
subspaces  in  a  p+m  dimensional  complex  vector  space  of  values  of  inputs  and  outputs.  Moreover, 
a  number  s  for  which  these  subspaces  intersect  nontrivially;  i.e.  dim(GsnKs)>0  is  a  closed  loop 
pole.  Following  this  line  of  thinking  Brockett  and  Byrnes  [21]  describe  a  Nyquist  stability 
criterion  for  the  general  p*m  case  based  on  frequency  response  data.  Their  result — which  is 
somewhat  abstract — provides  an  encirclement  condition  on  an  abstract  space  called  the  Grassman 
manifold.  Despite  the  abstract  nature  of  the  resulting  stability  test  it  does  retain  several  essential 
features  of  the  popular  single  loop  test  developed  by  Nyquist.  In  particular,  it  identifies  a  Nyquist 
contour  Fq  for  the  plant  as  the  image  of  the  imaginary  axis  s=jco  under  the  map  ker[G(s),-Ip]  and  a 
separate  object,  T^,  obtained  from  the  imaginary  frequency  response  of  ker[Im,K(s)]. 


In  [22]  Bennett  and  Baras  describe  a  geometric  stability  margin  based  on  an  angle  measure  of  the 
distance  between  the  intersection  of  the  respective  Nyquist  contours,  Tq  and  Tr.  Their  measure  is 
described  in  geometric  terms  but  is  meant  to  extend  to  the  general  MIMO  feedback  case  the  natural 
practical  notion  of  stability  margin  for  the  SISO  case  as; 

minks)+ f?y 

s=jto  I  1  Is/  I 

The  simplest  way  to  motivate  the  definition  of  the  geometric  stability  margin  is  to  recall  the 
definition  of  the  canonical  or  principal  angles  between  a  pair  of  subspaces. 

Definition  [24]:  The  principal  angles  0ke  [0^]  between  a  pair  of  complex  subspaces.X, 

Y,  under  the  assumptions,  dim  X  =p,  dim  Y=m  and  p<m,  are  given  recursively  for 
k=l,2,...,p  as, 

cose.  =  max  maxu'v  =  u*v 

K  ueX  veY  K  K 

ML*  H=1 

2  2 

subject  to  the  constraints 
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u*  u  =  0,  v*v  =  0 
J  J 

for  j=l,...,k-l.  Then  the  principal  vectors  for  the  pair  of  subspaces  are 

Canonical  angle  analysis  has  found  application  in  a  number  of  areas  including  the  computation  of 
statistical  correlations.  Numerical  algorithms  for  the  computation  of  canonical  angles  are  largely 
based  on  singular  value  decompostion  [24].  In  [24]  Bjorck  and  Golub  show  that  the  computations 
can  be  implemented  by  obtaining  a  pair  of  matrices  Q^»  Qy  of  dimension  nXp,  nXm, 

respectively  such  that  Q*  =  *  p ,  q,Qr  =  'm  and  Ox »  Qy  are  a  basis  for  X  (resp.  Y)— 
using,  for  example,  Gram-Schmidt  procedure.  Then  the  SVD  of  the  product; 

Q'xQY  = 

has  singular  values  as 

I  =  diag{a1,  ..,op}=cos0 

which  if  ordered  as  CT-j  -  •  •  •  -  ap  then  the  principal  angles  are  obtained  from 

©  =  diag{01,...>0p} 

with  V-*0P. 

Finally,  in  [22]  Bennett  and  Baras  show  that  via  an  alternate  form  of  the  above  principal  angle 
computation  one  can  associate  a  certain  minimum  singular  value  with  a  measure  of  how  nearly  two 
subspaces  intersect  in  a  nontrivial  way.  In  this  case  we  obtain  an  nX(n-p)  matrix  with 
orthonormal  columns  whose  span  is  the  orthogonal  complement  of  X.  Then  the  SVD  of  the 
product, 

q,  =  'vx  a  yy 

has  singular  values  related  to  the  canonical  angles  as 

A  =  sine. 
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Thus  the  minimum  singular  value  0  £  £  1  is  nonnegative,  real  number  which  indicates  how 

nearly  the  pair  of  subspaces  intersect.  Let  Y  =  X,(X,Y  )  ,  which  we  refer  to  as  the  minimum  gap. 


Definition:  The  geometric  stability  margin  for  the  closed  loop  system  is  given  as 
esm=ml"Y(£s.Es) 

S=j<0 


(2.2) 


In  [23]  we  show  that  the  computational  problem  of  orthonormalization  can  be  replaced  by  the 
introduction  of  a  pair  of  normalized  coprime  factors  for  G(s)  and  K(s). 

Definition:  A  pxm  transfer  function  G(s)  has  a  stable  normalized  right  (resp.  left)  coprime 
factorization  G=N  M'1  (resp.  G=lvHft)  if  N,M  (resp.  M,ft)  are  stable  transfer  functions 
which  satisfy, 

N*M+M*M=Im  (resp.  ft*ft+lM*lVi)=Ip).  (2.3) 

Moreover,  it  is  shown  in  [25]  that  every  pxm  transfer  function  has  a  Normalized  Right  (resp.  Left) 
Coprime  Factorization  (NRCF  or  NLCF)  which  is  unique  up  to  a  unitary  change  of  basis  in  input 
space  C  (resp.  output  space  ) . 

We  can  then  show  [23]  that  the  geometric  stability  margin  for  the  general  MIMO  case  can  be 
computed  as  follows.  Obtain  the  NLCF  of  the  plant  G=ftHR  and  the  NRCF  of  the  compensator 

K  =  PQ  .  Then  the  geometric  stability  margin  can  be  given  as 

9sm=m!n 
s=j«a 


(2.4) 


where 


=  ftp  +  MQ  (2.5) 

To  see  the  significance  of  this  construction  consider  that  the  normalization  processes  obtains  a 
frequency  dependent  change  of  basis  in  the  p-dimensional  space  of  outputs  such  that  the  matrices 


has  orth normal  columns  and 
[G,Ip]  =SH[Rffl] 
has  orthonormal  rows. 
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An  alternative  expression  for  plant  uncertainty  advocated  by  Vidyasagar  [25-27]  in  terms  of 
additive  stable  perturbations  to  the  factors  in  a  coprime  factorization  of  the  plant  has  certain 
advantages  in  analysis  of  feedback  systems.  In  recent  work  Glover  and  McFarlane  [31]  also 
considered  this  class  of  unstructured  perturbations  and  obtained  a  surprisingly  explicit  and 
intuitively  appealing  solution  to  the  corresponding  robust  stabilization  problem  when  the  coprime 
factorization  is  normalized.  In  particular,  they  considered  the  problem  of  stabilization  of  a  nominal 
plant  transfer  function  G(s)  and  a  family  of  transfer  functions  given  in  terms  of  its  NLCP 


G^-1^ 


given  as 


Ge=  { (M+Am)'  1  (54+An):  ll[ltfR]1loo<e} 


where  II. IL,  is  the  H°°  norm  of  the  transfer  function  given  as 


(2.6) 


||G^)|!oo=supama)j[G(s)} 

s=jco 

and  Am  ,  An  are  stable  unknown,  transfer  functions  which  represent  the  uncertainty  in  the 
nominal  plant  model.  It  was  demonstrated  by  Vidyasagar  [25-27]  that  this  description  of  plant 
uncertainty  has  advantages  over  additive  or  multiplicative  unstructured  uncertainty  models.  For 
example,  the  number  of  unstable  poles  may  change  as  the  plant  is  perturbed.  Figure  2  illustrates 
the  class  of  systems  described  by  perturbations  of  normalized  coprime  factors. 


Figure  2:  Perturbation  of  NCPF 


It  is  now  apparent  the  role  of  coprime  factorization  originally  conceived  in  the  work  of  Youla, 
Rosenbrock,  and  others  was  somewhat  constraining.  The  goal  of  much  of  this  reaserch  was  to 
extend  the  obvious  constructions  for  rational  scalar  transfer  functions  of  factorization  in  terms  of 
numerator  and  denominator  polynomials  to  the  required  matrix  constructions  for  multiloop 
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systems.  More  recently,  connections  with  state  space  constructions  have  become  evident  by 
introducing  an  alternate  viewpoint  [25].  Instead,  we  think  of  coprime  factorization  over  an 
alternate  algebraic  ring  of  rational  tranfer  functions  whose  poles  are  contained  in  some  arbitrary  left 
half  plane  within  C.  Since  in  most  cases  we  are  concerned  with  asymptotic  stability,  it  suffices  to 
consider  coprime  factorizations  over  the  space  of  all  rational  functions  of  the  Laplace 

variable  s  analytic  for  Re  s>0.  Then  a  state-space  construction  for  the  normalized  left  (resp.  right) 
coprime  factorizations  over  51 can  be  obtained  in  terms  of  the  solution  to  the  Control  (resp. 


Filter)  Algebraic  Riccati  Equation: 


ATX+  XA  -  xbbtx+  ctc  =  0 
AY  +  YAT  -  YCTCY  +  BBT  =  0 


(CARE) 

(FARE) 


Specifically,  let  G(s)  be  realized  by  a  linear  state  space  system; 


*(t)  =  A$(t)  +  Bu(t) 

z(t)  =  cm 

so  that  Q(s)  =  c[s  I  -  A]  B  with  (A,B,C)  a  minimal  realization.  Then  the  respective  left 
normalized  coprime  factors  can  be  shown  to  be  realized  as 


fi(s)  =  C[sI-A0]_1B 
S4(  s)  =  I-CfsI-Aol^H 


(2.7a) 

(2.7b) 


where 


H=  -  YC  ,  A  =  A  +  HC 


The  geometric  stability  margin  provides  a  numerically  attractive  measure  of  closed  loop  stability 
which  includes  the  class  of  perturbations  just  described.  It  will  recognize  possible  pole/zero 
cancellations  which  might  occur  in  the  formation  of  any  variation  of  the  return  difference  matrices 
and  it  provides  a  characterization  of  plant  neighborhoods  for  robust  analysis  as  above  in  terms  of 
metric  specifically  designed  for  the  stability  margin  for  multiloop  stability  directly  in  terms  of  the 
generalized  Nyquist  criterion  [21]. 

2.2.  Design  Methods  for  Optimal  MIMO  Stability  Margins 

The  primary  reason  for  the  introduction  of  coprime  factorization  by  Youla  et  al  [28]  was  to  facilitate 
analytic  optimization  for  the  resolution  of  engineering  tradeoffs  in  design.  The  coprime 
factorization  permits  optimization  to  be  carried  out  without  regard  for  stability  of  the  closed  loop 
system.  This  fact  follows  from  the  following  well  known  lemma. 
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Lemma:  For  each  real  rational  G(s)  there  exists  left  and  right  coprime  factorizations; 
G=NM_1=  each  contained  in  and  functions  V.U.V.OeSRH^  such  that 


'V  -Tj  Y  M  U  - 

-10- 

.  -Tift  J-  N  V  . 

.01. 

With  the  above  notation  the  Youla  parametrization  provides  a  complete  parametric  description  of  all 
stabilizing  compensators  for  a  given  plant  transfer  function  G(s). 

Lemma:  Any  compensator  transfer  function  K(s)  which  stabilizes  the  plant  G(s)  in  closed 
loop  can  be  factored  as 

K  =  (U+MQ)(V+NQ)-1  =  (V+SQ)-l(0+ltfQ)  (2. 10) 

for  some  Qe 

This  has  been  called  the  Q-parametrization  by  Desoer  et  al  [60].  It  is  significant  to  recognize  that 
this  parametrization  guarantees  that  the  closed  loop  system  (see  Figure  2)  is  internally  stable-,  i.e  . 
the  composite  transfer  function  satisfies, 

K(,p-Vl  h_ 

G  (lm-KG)  . 

The  focus  of  the  current  phase  1  project  is  to  investigate  the  application  of  some  recent  results  of 
K.  Glover  and  D.  McFarlane  [31-33]  on  robust  stabilization.  In  the  next  few  paragraphs  we 
summarize  the  basis  for  their  results. 

Defintipn:  A  feedback  system  consisting  of  plant  transfer  function  G(s)  and  compensator 
KCs)  is  robustly  stable  with  respect  to  the  margin  e  (or  e- robustly  stable)  if  and  only  if  the 
feedback  system  of  (2. 1)  K,  Ge  is  internally  stable  for  any  Ge  e  Ge  (as  given  in  (2.6)). 

The  symbol  e  represents  a  stability  margin  for  the  system  in  the  sense  that  it  indicates  the  'size'  of 
tolerable  perturbations  internal  to  the  plant  (in  terms  of  the  coprime  factorization).  (As  we  will 
demonstrate  shortly,  this  stability  margin  is  identical  to  the  geometric  stability  margin  deser  !x\ 
above.)  The  robust  design  problem  is  to  stabilize  not  only  the  nominal  plant,  G,  but  the  family  of 
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perturbed  plants  defined  in  (2.6)  using  a  single  feedback  controller  K(s)  (see  Figure).  Their  first 
result  provides  necessary  and  sufficient  conditions  for  robust  stabilization  of  a  given  plant. 

Lemma:  The  feedback  system  with  plant  G=lvH5J  is  e-robustly  stabilizable  by  some 
compensator  K  if  and  only  if 


inf  j|  K(LGKy1M1  |]  l 

K(s)e9moc  (i-GKyW  .  ~  ~e 


(2.11) 


This  is  stated  in  the  form  of  an  H°°  optimization  problem  which  could  be  solved  by  the  standard 
iterative  procedures  outlined  in  Francis  [20].  To  follow  this  approach  one  could  first  proceed  by 
recasting  the  optimization  problem  as  a  “standard  problem”.  This  is  illustrated  in  Figure  3. 


Figure  3:  Robust  Stabilization  Problem  in  Standard  Form 


This  can  be  seen  by  comparing  the  standard  problem  of  Figure  1  where  in  this  case  we  take 


P  =  |_M  J  IVN 
M'1  M_1N 


LM  M  N  j,  (2.12) 

Then  the  standard  H°°  problem  is  to  minimize  the  H°°  norm  of  the  transfer  function  from  w  to  e; 


(2.13) 


P11+P12K«-P2^>' 


2  1 


K(I-  GK)_1  M"f 
.  (I  -  GK)_1M_1. 


which  demonstrates  the  interpretation  displayed  in  Figure  3. 


However,  by  using  normalized  coprime  factorizations,  Glover  and  McFarlane  show  that  the 
iterative  procedures  can  be  avoided  altogether,  and  an  exact  solution  can  be  obtained  very  simply. 
The  key  idea  is  to  combine  the  Youla  parameterization  of  stabilizing  controllers  described  above 
together  with  interpolation  results  of  Nehari  and  Glover  for  the  computation  of  the  Hankel  norm. 
In  particular,  they  show  that  the  plant  G  is  e-robustly  stabilizable  if  and  only  if 


e2<l-lllM,ft]H„  (2.14) 

where  the  Hankel  norm  of  a  stable,  proper  system  G(s)=C[sI-A]_1B  is  given  as 

l|S^)lH=  =  .  (2  15) 

i.e.,  the  maximum  singular  value  of  the  associated  Hankel  operator, 


(r^(t)  =  /CeA<f+T)Bv(t)dt 

o  .  (2.16) 

It  is  also  shown  in  [30]  that  the  Hankel  norm  can  be  obtained  from  the  state  space  realization  of  the 
stable  system  via  the  solutions  of  the  pair  of  associated  Lyapunov  equations; 

AP+  PA*  +  BB*  =  0 

AtQ+QA  +  C,C=0,  (2.17) 

as  shown  above. 


Using  these  constructions  Glover  and  McFarlane  show  that  the  optimal  stability  margin  for 
G(s)=C[sI-A]1B  given  by  the  minimum  achievable  e0  in  (2.1 1)  can  be  directly  obtained  from  X, 

Z>0,  the  positive  definite  solutions  of  the  pair  of  Riccati  equations  (CARE)  and  (FARE)  as 

°  yi+x-maj(ZX)  .  (2. 1 8) 

Significantly,  Glover  and  McFarlane  also  provide  formulae  for  the  characterization  of  all 
controllers  which  achieve  e-robust  stabilization  and  from  the  previous  discussion  such  level  is 
achievable  only  if  e<eo,  with  eq  the  optimal  stability  margin  obtained  from  (2.14).  Formulae  for 
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the  desired  e-robust  compensators  are  given  directly  in  terms  of  state  space  realizations  which 
facilitates  computer  based  computation.  Details  are  given  in  [32].  It  is  also  potentially  important 
to  recognize  the  flexibility  available  in  this  characterization  which  permits  consideration  for  various 
additional  design  tradeoffs.  In  this  study  we  have  focused  on  the  simplest  form  for  realizing  an  e- 
robust  compensator.  This  is  summarized  as  follows. 


Theorem:  A  compensator  K(s)  which  achieves  the  specified  stability  margin  e,  0<£<£o,  is 
obtained  in  the  form 


KIsHC^sl-AjX 


(2.19) 


where  the  state  space  realization  is  given  by  the  matrices 

Ak= Ac+Y2(wV1zc,a 
Bk  =  Y^wV’zc', 


and  where  the  matrix  terms  are 

Ac  =  A  +  BF, 


w=  xz-g, 


(2.21) 


and  the  scalar  terms  are _ 

y~Tt  C  =  v Y2  “  V  (2.22) 

The  following  demonstrates  the  connection  between  the  geometric  stability  margin  and  the  e 
stability  margin  used  above. 


Claim:  The  geometric  stability  margin  is  identical  with  the  e  stability  margin  introduced  by 
Glover  and  McFarlane  in  the  context  of  the  robust  stabilization  problem. 

Proof:  To  see  this  consider  a  plant  transfer  function  together  with  NLCP  factorization  G=  M 
over  Si  H  °°.  Then  for  some  K(s)  which  stabilizes  G(s)  the  achieved  stability  margin  e  is  defined 

in  terms  of  the  H°°  norm  as, 
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Y  =  t  = 


K 

L*  PJ 


(lp-GK)"V 


=  SUP  oma] 

s=jto 


TK_ 

J  P. 


-1  --0 

(lp-CK)  M  (s) 


Let  K  -  PQ  be  NRCP  factorization  of  the  compensator  over  H 80 .  Then  it  is  clear  that 


K 

l'  p. 


(Ip  -  GK)  'm  1  = 


P 

LQI 


[MQ-ra3] 


and  since  is  an  inner  function  the  resulting  norm  can  be  obtained  as 

J.=  I[MQ-hP]_l1 

or  equivalently, 

£  = 1  n  * 

S=J<U 

where 

<D(s)  =  [MQ-NP](s) . 

Thus  we  conclude  that  e  =  9Sm- 


The  significance  of  the  geometric  stability  margin  is  to  lend  an  insight  relating  the  classical  Nyquist 
analysis  to  the  MIMO  setting  and  to  clarify  the  role  of  the  normalization  in  geometric  terms.  Thus 
it  is  clear  from  the  definition  of  the  geometric  stability  margin  that  0sm  is  the  sine  of  the  minimum 
principal  angle  between  a  certain  pair  of  subspaces  and  therefore  0<e<l.  We  suspect  that  the  role 
of  the  geometric  picture  may  offer  new  insights  in  terms  of  frequency  shaping  design  for  combined 
stability  and  performance  robustness.  In  the  Phase  1  study  we  have  highlighted  these  features 
briefly  as  part  of  the  flight  control  benchmark. 

2.3.  Design  Methods  for  combined  performance  and  stability  robustness 

Loop  shaping  is  a  well  established  design  methodology  for  control  design  for  MIMO  systems  with 
LQG/LTR  being  one  of  the  better  developed  techniques.  An  advantage  of  loop  shaping  is  that  a 
suitable  loop  shape  for  both  stability  and  performance  robustness  can  often  be  chosen— at  least  in 
the  case  of  unstructured  model  uncertainty— with  relative  ease.  The  motivation  for  loop  shaping 
comes  from  classical  frequency  domain  design  for  SISO  systems  where  lead/lag  circuits  as  well  as 
PID  circuits  provide  the  component  parts  for  attaining  the  desired  loop  shapes  and  stability  margins 
guide  the  designer  in  quantitative  choices.  For  example,  following  SISO  design  we  expect  that 
low  frequency  gains  should  be  sufficiently  high  for  typical  performance  requirements  of  good 
tracking  at  steady  state,  while  high  frequency  gain  should  be  rapidly  attenuated  in  regions  where 
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the  model  dynamics  may  neglect  significant  parasitics.  The  achievable  roll  off  in  SISO  systems  is 
related  to  phase  in  a  direct  way  as  described  by  Bode.  For  MIMO  systems  it  is  in  general  also 
difficult  to  guarantee  good  robustness  in  the  gain  cross  over  region — particularly  for  nonminimum 
phase  plants. 

A  new  method  of  loop  shaping  has  been  proposed  by  McFarlane  and  Glover  [33]  and  has  been 
examined  for  the  flight  control  design  considered  in  this  study.  In  this  approach,  the  results  of  the 
problem  of  optimal  robust  stabilization  of  the  normalized  coprime  factors  (i.e.,  maximization  of  the 
geometric  stability  margin)  are  used  as  a  basis  for  design.  As  an  example  consider  the  robust 
stabilization  of  a  simple  SISO  plant  and  choose  arbitrary  stable  NLCP  factors; 

Qs)  =  l  M(s)  =  ^.  Nts)  =  jV 

and  take  as  compensator 

K(s)  =  -  k . 

The  class  of  systems  that  can  be  stabilized  by  K(s)  can  be  described  as 


g.<s>  =  (m+am)'  (n+an) 

with  ^  where  ^  is  the  geometric  stability  margin  which  can  be 

computed  as 

0sm(G,K)=1/ 


K1 


(l-GK)  V1 


In  this  case  we  can  obtain 


0sm(l/S,  -k)  = 


(1+  k  2)  1/2,  for  k<1 
(1+  k2)  1/2,  f  o  r  k  >  1 


Hence  the  stability  margin  is  poor  if  either  k«l  or  k»l  and,  in  fact,  optimal  if  k=l.  In  the 
former  case  the  closed  loop  pole  is  moved  to  -k  and  hence  a  small  change  in  the  open  loop  pole — 
say  to  +2k — would  not  be  stabilized  by  the  fixed  feedback.  In  the  latter  case,  the  gain  is  high  so 
that  phase  errors  could  cause  loss  of  stability.  Note  that  in  both  cases  the  interpretation  just 
described  will  depend  on  several  factors,  including,  for  example,  the  scaling  of  time,  inputs,  and 
outputs.  Note  also  that  the  loop  shape,  I/O),  has  been  substantially  altered  in  either  case  if  k«l  or 
k»l. 


The  method  proposed  in  [33]  and  implemented  in  this  study  is  as  follows: 

(i)  Loop  shaping :  Choose  weighting  functions  Wj  and  Wo  such  that  the  weighted  or 
shaped  p!antGs=WoGWj  has  the  desired  loop  shape. 

(ii)  Robust  Stabilization :  Compute  the  optimal  stability  margin  eo=m^xBsm(Gs,K) 

as  in  (2.18).  If  eo(Gs)«l,  return  to  (i),  otherwise  compute  a  controller  Ks  which 
achieves  the  maximal  stability  margin  for  the  weighted  plant  using  (2. 19) — (2.22). 

(iii)  Control  Realization :  The  final  controller  is  realized  by  K=WjKsWo. 

The  design  method  proposed  has  several  useful  and  assured  properties.  For  example,  the  direct 
optimization  of  the  stability  margin  assures  the  geometric  stability  margin  is  achieved  or  execeeded 
at  all  frequencies.  This  feature  is  demonstrated  in  the  examples  considered  in  this  report  and 

provides  a  graphic  demonstration  of  the  broadband  matching  properties  of  H°°  optimization. 
Second,  in  the  high  gain  region  the  controller  Ks  computed  for  the  shaped  plant  reduces  the  loop 
gains  (i.e.  min/max  singular  values)  by  a  factor  of  at  most  e0  (^s)  .  Third,  in  the  low  gain  region 
the  controller  Ks  can  increase  the  loop  gains  by  at  most  a  factor  of  l/eo(Gs)  Moreover,  if 
eo(Gs)«l,  then  0sm^eo(Gs)«l,  for  any  stabilizing  controller.  Hence,  in  this  case  the  loop  shape 
is  not  compatible  with  the  robust  stabilization  requirement  and  any  controller  which  achieves  n  will 
have  a  poor  stability  margin  in  the  gain  cross  over  region,  or  undesirable  decrease  in  the  low 
frequency  gain  (i.e.,  loos  of  perfromance),  or  undesirable  increase  in  the  high  frequency  gain  (i.e., 
loss  of  stability  margin). 

We  therefore  interpret  the  optimal  stability  margin  as  an  indicator  or  engineering  figure  of  men;  tor 
the  compatibility  of  the  loop  shaping  proposed  in  step  (i).  The  resulting  performance  will  be 
robust  to  any  small  changes  in  plant  parameters  which  can  be  absorbed  into  the  unstructured  model 
uncertainty  in  terms  of  the  perturbations  to  the  NLCP  factors.  However,  large  changes  in  real 
parameters  arising  from  changes  in  the  system  operating  point  may  be  difficult  if  not  impossible  to 
embed  in  perturbations  to  the  NLCP  factors  of  the  plant  without  making  the  design  overly 
conservative.  An  example  of  this  fact  is  discussed  in  the  flight  control  benchmark  considered  in  a 
later  section. 

3.  Worst  Case  Control  Design  for  Combined  Structured-Unstructured 
Model  Uncertainty 


In  Figure  3. 1  we  display  an  uncertain  system  with  feedback  controller,  K,  A(s)  is  the  plant  model 

uncertainty  and  it  is  desired  that  the  transfer  function  from  w  to  z  remain  bounded  less  than  1  for  all 
A  such  that  11^  1  .  The  bound  can  be  rescaled  to  any  value  and  this  problem  can  include 

both  performance  and  stability  robustness. 


A  sufficient  condition  for  robust  performance  is  that  the  closed-loop  transfer  function  from 

w  =  [v  ,  W  ]  to  2  =  [e  i  z  ]  have  an  H°°  norm  less  than  1.  This  condition  is  not  conservative 
if  robust  stability  or  nominal  performance  alone  are  considered,  but  will  become  conservative 


when  robust  performance  is  required  [16].  H°°  optimization  is  an  effective  tool  which  can  be 

applied  to  such  problems  to  maximize  the  robustness.  Additionally,  the  computational 
requirements  are  now  understood  [35].  The  resulting  approach  is  computationally  quite  tractable 
and  provides  controllers  with  degree  no  greater  than  the  degree  of  the  model  P(s).  However,  it 
does  not  address  structured  uncertainty. 


Figure  3.1:  Uncertain  Structure  of  Feedback  Control  Problem. 


3.1.  Summary  of  Some  Existing  Methods  for  Design 
3.1.1.  Structured  Singular  Values 

The  robust  performance  problem  of  Fig.  3.1  for  1  and  subject  to  a  specific  block 

diagonal  structure  for  A  is  considered  by  Doyle  in  terms  of  the  'structured  singular  value'  analysis. 
For  a  large  number  of  blocks  this  problem  is  not  tractable,  but  a  good  approximation  can  be 
obtained  by  introducing  diagonal  'D-scales'  on  the  inputs  v  and  outputs  z.  These  D-scales  are 
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chosen  to  be  functions  of  frequency  that  minimize  the  H°°  norm  of  the  scaled  transfer  function 
W  ->?.  This  can  provide  an  effective  analysis  tool. 

The  design  problem  is  more  difficult.  It  can  be  addressed  by  alternately  minimizing  over  the  D- 
scales  and  then  over  the  stabilizing  controllers  K  until  convergence.  There  are  few  theoretical 
results  to  justify  this  procedure.  Real  parametric  variation  can  (in  principle)  aLo  be  addressed 
using  a  nontrivial  modification  of  the  method  [36].  Again,  this  is  an  analysis  tool  and 
computationally  more  involved  than  the  case  of  complex  (i.e.  frequency  dependent)  structured 
uncertainties. 

3.1.2.  Methods  Based  on  Kharitonov's  Theorem 

Kharitonov's  theorem  provides  a  remarkably  efficient  test  for  the  stability  of  a  class  of  systems 
when  the  closed  loop  characteristic  polynomial  can  be  written  as  a  polynomial  with  coefficients— 
each  of  which  lies  in  an  independent  interval  of  the  parameter  space.  This  result  has  lead  to  a 
variety  of  methods  for  analysis  of  interval  matrices  whose  elements  are  parameters  contained  in 
known  intervals.  However,  for  more  general  structures  the  resulting  tests  become  intractable 
becouse  of  the  size  of  the  test  matrices  involved.  Such  tests  are  based  on  determinantal  inequalities 
and  involve  purely  algebraic  computations  whose  complexity  leads  to  algorithms  which  are  NP- 
complete  (in  the  jargon  of  computer  science.) 

3.1.3.  Quadratic  Stabilization 

The  work  of  Peterson  and  Hollot  [37]  on  quadratic  stabilization  of  uncertain  systems  considers 
norm-bounded  but  possibly  time-varying  perturbations  of  a  nominal  system.  Furthermore,  they 
employ  a  strong  stability  condition  which  requires  the  existence  of  a  quadratic  Lyapunov  function. 
It  can  be  shown  that  this  version  of  stability  is  identical  to  requiring  an  H°°  norm  bound  on  the 
transfer  function  from  the  outputs  of  the  perturbations  back  to  the  corresponding  inputs  (as  in  Fig 
3.1).  Robust  stability  then  follows  from  a  small  gain  argument.  Results  using  this  procedure  can 
therefore  be  interpreted  in  an  H°°  robust  stabilization/performance  framework. 

3.1.4.  Combined  LQG-H°°  Design 

Bernstein  and  Haddad  [38]  consider  a  combined  LQG  and  H°°  control  design.  Their  introduction 
suggest  that  they  are  minimizing  the  Hz  norm  subject  to  a  constraint  on  the  H°° norm.  This  is  in 
fact  not  the  case  as  a  careful  examination  of  the  paper  reveals.  Instead  they  minimize  upperbounds 
rather  than  the  actual  norms  as  claimed.  In  many  circumstances  this  will  result  in  a  good 

approximation  to  the  original  problem.  Indeed,  in  the  case  when  the  H°°  constraint  and  the  Hz 
criterion  are  both  applied  to  the  same  closed  loop  transfer  function,  then  (as  shown  by  Mustafa 
( 39])  their  criterion  is  identical  to  the  entropy  maximization  criterion  in  H°°  design  as  considered  by 
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Mustafa  and  Glover  [40]  (see  also  appendix).  The  entropy  criterion  for  a  closed  loop  transfer 
function,  H(s),  is 

oo 

I(H,y)  =  -— 2  1  log[det{  l-y2H(jo))H*(jco)}]do)  (3.1) 

*  -OO 

Clearly,  it  is  required  that  IIHIL^y  for  I(H,y)  to  be  well-defined  and  it  is  easily  shown  that  "\/ I(H,y) 

2 

is  an  upper  bound  on  the  H  norm. 

The  major  contribution  of  these  papers  is  to  allow  different  transfer  functions  in  the  and  H°° 
criteria  and  can  therefore  address  the  problem  of  reduced  order  compensator  design,  the  main 
disadvantage  being  the  substantial  computational  difficulty  involved  in  solving  the  resulting 
coupled  Riccati  equations. 

In  [41]  Yeh  et  al  consider  the  problem  of  designing  controllers  for  the  combination  of  both  real 
parametric  uncertainty  and  unmodeled  dynamics.  The  approach  utilizes  a  combination  of  the 
methods  of  Bernstein  and  Haddad  together  with  the  results  of  Peterson  and  Hollot.  However, 
since  both  approachs  employ  upper  bounds  one  can  expect  the  results  to  be  conservative. 

To  be  specific  consider  the  closed  loop  system 

x  =  (A  +  AA)x  +  Bw 
z  =  Cx 

where 

P 

AA=£DjAjEj,  O^Sl. 

i-1 

This  system  can  be  rewritten  in  the  form  of  Figure  3.1  by  defining, 

e.  =  Ex,  v.=A.e. 

i  i  i  ii 

to  give, 

P 

x  =  Ax  +  £D.  v.  +  Bw 

i=l  '  1 

The  corresponding  Riccati  equation  employed  in  [41]  is  then 

0  =  AQ  +  QA*  +  y^Qc'cQ  +  Q  +  BBl 


(3.2) 

(3.3) 

(3.4) 

(3.5) 
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where 


P  (  P 

n  -  £D.  d  +  C  IE.  E1. 

i=1  Vi  =1 


Q 


(3.7) 


However,  the  existence  of  Q>0  to  satisfy  (3.7)  is  precisely  a  test  on  the  H°°  norm  ot  the  transfer 
function  from  w  to  z  in  Fig.  3.1.  Hence  the  performance  will  be  robust  to  the  uncertainty  as 
claimed  but  will  also  be  robust  to  the  unstructured  uncertainty. 

This  approach  gives  a  viable  design  technique  for  structured  and  unstructured  uncertainty, 
however,  it  will  be  conservative  because  it  essentially  (although  not  obviously)  embeds  the 
structured  uncertainty  into  a  larger  unstructured  uncertainty  perturbation  (e.g.  oiock  ol’ A),  in  iaci 
if  the  coprime  factors  Mj  and  Nj  are  optimally  scaled  with  frequency,  then  thL  scheme  become* 

identical  to  the  structured  singular  values. 

3.2.  Worst  Case  Design  as  Minimax  Optimal  Control  Problem 

The  geometric  stability  margin  0sm(Gs,K)  has  been  shown  to  be  a  suitable  design  indicator  and 

engineering  figure-of-merit  for  both  performance  and  stability  robustness  of  closed  loop  control 
systems.  However,  the  control  design  procedure  suggested  above  will  be  conservative  if 
structured  uncertainty  due  to  real  parametric  variation  is  present  in  the  model.  Instead,  in  this 
section,  we  take  the  plant  model  to  be  of  the  form  Ga(s)  =  G(S,a)  ,  a  function  of  a  real  parameter 

k 

vector  a.  If  we  assume  a  e  0  s  9?  a  compact  set  in  the  k-dimensional  real  vector  space  then 

the  optimal  ’’robust"  design  solution  for  stability  is  to  find  a  compensator  K  which  staoilmos  ihe 
nominal  plant  GQ(s)  and  maximizes  the  worst  case  stability  margin  attained  over  the  >e: 
k 

a  e  0  q  Jt  j-fere  let  y(a,K)=l/0srn(Ga,K)  then  we  wish  to  find  K  which  solves; 

inf  sup  Y(a,  K) 

K  ae©  .  (3.8) 

Io  see  formally  that  this  is  appropriate  for  combined  structured  and  unstructured  uncertainty 
consider  the  problem  as  follows:  given  the  design  model 

(plant)  y  =  G(a)u 
(comp)  u  =  Ky 


with  u  an  m-vector  of  plant  inputs  and  y  a  p-vector  of  plant  outputs.  See  Figure  3.2.  The 
structured  plant  uncertainty  is  characterized  by  a  I-vector  of  parameters,  ae  0,  with  20. 


Figure  3.2:  Robust  Control  Design  Model  for  Mixed  Uncertainty 

Problem  of  Robust  Stabilization:  Find  K(s)  mXp  (controller)  fixed  wrt.  a  which  stabilizes  a  class 
of  plants  G(a)  given  by: 

(uc)  Unstructured  uncertainty  conditions:  for  each  ae  0,  G  =  is  taken  to  be  a 
normalized  left  coprime  factorization  (NLCP);  i.e.,  (N,M)  are  left  coprime  and 
MM*  +  NN*  =  I 

G€Ge  =  {  (M+  Am)_1  (N  +  An)  :  ||  [Am  ,  An]  ||  _<  e)  (3.9) 

(sc)  structured  uncertainty  conditions: 

GeGa  =  {  G(s;  a) :  oce0  },  ft1  2©.  (3.10) 

For  the  mixed  (or  combined)  uncertainty  conditions  we  simply  mean 
Ge  G£  n  Ga  . 


A  closely  related  but  more  tractable  problem  is  the  maximin  (or  dual)  version  of  (3.8); 


sup  inf  Y(a,K)  = 

ae  0  K 


(3.11) 


whose  solution  is  now  feasible  because  of  the  results  of  Glover  and  McFarlane  [32]  as  described  in 
Section  2.2  of  this  report.  These  results  provide  a  simple  and  computationally  straightforward 
solution  of  the  first  level  optimization  problem;  i.e.,  the  optimal  stability  margin  can  be  computed 
for  any  choice  of  a.  The  use  of  minimax  design  for  parametric  model  uncertainty  has  been 
considered  previously  in  only  a  few  cases  which  have  been  reported  in  the  open  literature. 
Examples  available  provide  an  insufficient  basis  for  theoretical  extension  of  the  results  to  large, 
complex  problems  but  do  suggest  the  practicality  for  relatively  simple  systems  with  small  number 
of  parameters  [43]. 


Sufficient  conditions  for  the  minimax  and  maximin  problems  to  have  the  same  saddle  point 
solution  include  the  function  y(a,  1^  being  convex  in  K  and  concave  in  a  e  0.  It  is  not 
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apparent  that  either  of  these  conditions  are  likely  to  be  satisfied  globally.  It  can  however,  be 
shown  that  if  we  employ  the  Youla  parametrization  for  the  choice  of  K  (i.e.,  to  find  a  single 
controller  for  all  a  e  0 )  then  the  objective  will  be  convex  in  the  Q  parameter  [20],  The  behavior 
with  respect  to  oc  e  0  is  unlikely  to  be  concave  except  locally.  At  a  saddle  point  it  can  be  checked 
whether  ®  sm^a’  is  locally  convex  with  respect  to  a.  To  do  this  we  find  the  relation  KQ(a) 

by  the  minimization  in  (3.1 1);  i.e.,  the  optimal  controller  depending  on  the  unknown  parameters 

Then  at  the  point  (ctQ,  we  check  that  ®sm^a»  is  minimized  with  respect  to  a  e  0. 
Furthermore,  if  we  use  the  suggested  compensator  KQ  then  we  have  the  guaranteed  bound 


and  we  can  evaluate  the  achieved  performance 
suP8sm<Ga.Ko> 

as  0  (3.1.1) 

to  establish  the  sadd’“  point  equivalence  by  comparison  with  the  lower  bound  (3.12)  The 
combination  of  the  steps  (3.11)  and  (3.13)  will  in  general  require  extensive  evaluations  of  the 
stability  margins  at  various  candidates  (K,a)  which  may  create  extensive  computational  burden 
However,  as  illustrated  in  the  benchmark  example,  several  significant  flight  control  design  issues 
can  be  addressed  using  this  approach.  We  emphasize  also  that  since  no  bounding  approximations 
are  used,  the  structure  of  the  parametric  uncertainty  is  retained  in  an  essential  way  in  the  design 
process.  It  is  also  apparent  that  the  worst  case  design  obtained  can  be  directly  understood  from  the 
model  assumptions  and  the  objective  of  a  robust  solution  [42], 

The  minimax  and  maximin  problems  will  also  be  equivalent  in  the  case  when 

suPesm(Ga'K)  =e<Ga  -  K> 

ae  0  0  for  all  K;  (3.14) 

i.e.  when  the  worst  case  set  of  parameters  is  independent  of  the  choice  of  controller.  This  is  often 
the  case,  for  example,  when  the  extreme  values  of  the  parameter  correspond  to  maximum  phase  lag 
or  ultimate  instability  (see  example  in  the  next  section). 

4.  Control  Design  Benchmark-Robust  Control  for  Aircraft  with 
Relaxed  Static  Stability 


4.1.  Longitudinal  Dynamics 


We  summarize  the  basic  equations  of  motion  which  govern  the  longitudinal  dynamics  of  an 
aircraft.  Further  details  may  be  found  in  [49-52],  Figure  4.1  identifies  the  body  axes  (X-Z)1,  the 
velocity  V,  the  body  attitude  0,  the  flight  path  angle  y,  the  angle  of  attack  oc=0+y,  and  the  principle 
forces  acting  on  the  airframe.  These  include  lift,  drag,  thrust  and  weight. 


Figure  4.1.  Principle  Longitudinal  Parameters 

The  basic  equations  of  motion  include:  linear  momentum  balance  in  the  X  direction,  linear 
momentum  balance  in  the  Z  direction  and  angular  momentum  balance. 

m(u+w0)  =  -mgsin0  +  Lisina  +  Ltsinat  +  T  -  Dcosa  (4. 1  a) 

m(w  -u0)  =  mgcos0  -  Lwcosa  -  Ltcosoq  -  Dsina  (4. lb) 

10  =  Mw  +  lwLwCosa  -  ltLtCosoq  -  c0  (4.1c) 

where  at  is  the  tail  angle  of  attack  and  is  related  to  the  angle  of  attack  a,  pitch  rate  0,  tail  angle  it, 
downwash  angle  e  and  the  elevator  deflection  angle  5e  via  the  relation 

at  =  a  +  it  -  e  +  8e  +  Ot/V)0 

We  also  have 

lw  +  lt  =  l* 

'The  X  axis  is  usually  aligned  so  that  o=0  corresponds  to  /.cro  lift. 
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The  lift  and  drag  forces  depend  on  the  velocity  V,  air  density  p  and  surface  area  S  via  the  relations 


L  =  CL(a)^pV2S 

(4,?a) 

D  =  CD(a)^pV2S 

v4.2b) 

M  =  CM(ct)^pV2S 

1,4.2c) 

In  view  of  the  dependence  of  the  forces  and  moments  on  V,  a  it  is  convenient  to  replace  u,  w  in 
(4.1)  by  V,  a  using  the  transformation  relations 


w  =  Vsina  (4.3a) 

u  =  Vcosa  (4.3b) 


Thus,  we  obtain 

m(cosaV-Vsinaa  +Vsina0)  =  -mgsinG  +  Lwsina  +  Ltsinat  +  T  -  Dcosa  (4.4a) 

m(sinaV+Vcosaa  -Vcosa0)  =  mgcos0  -  Lwcosa  -  LtCosat  Dsina  (4.4b) 

10  =  Mw  +  lwLwCosa  -  ltLtcosat  -  c0  (4.4c) 


Notice  that  equations  (4.4)  can  be  organized  in  the  vector  form 


mcosa 

-m  Vsina 

mVsina 

0 

" 

"  V  " 

msina 

mVcosa 

-mVcosa 

0 

d_ 

a 

0 

0 

1 

0 

dt 

e 

L  0 

0 

0 

I 

- 

_  0  - 

-mgsin0  +  Lwsina  +  Ltsinat  +  T  -  Dcosa  *”| 
mgcos0  -  Lwcosa  -  Ltcosat  -  Dsina 
0 

Mw  +  lwLwcosa  -  ltLtcosat  -  c0  _J 


(4.5) 


4.1.1  Nondimensional  Equations  of  Motion 
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Let  us  introduce  a  normalized  velocity  by  identifying  a  nominal  velocity  (for  example,  the 
maximum  cruise  velocity)  V0  and  define 

V  :=  (V/V o) 


and  also  the  nondimensional  quantities 
k  :=  (lw/1*) 

Aw  :=  (Lw/mg),  At  :=  (Lt/mg),  A  :=  (D/mg),  fl  :=  (T/mg),  Sw  :=  (Mw/l*mg) 


Equation  (4.5)  can  now  be  written 


cosa 

-■osina 

osina 

0 

— 

o 

sina 

ocosa 

-ocosa 

0 

d 

a 

0 

0 

1 

0 

dt 

0 

L  0 

0 

0 

1 

- 

-  (Vo/g)0  - 

-sin0  +  Awsina  +  Atsinat  +  n  -  Acosa 
cos0  -  Awcosa  -  Atcosat  -  Asina 
(Vo/g)0 

(VQl*/gr2){5;w  +  KAwcosa  -  (l-K)AtCosat)  -  (cVo/mgr2)(vo/g)0 


Now,  let  us  introduce  a  nondimensional  time,  x,  and  pitch  rate,  q, 
T  :=  (g/Vo)t,  q  :=  (Vo/g)0 


in  order  to  obtain  the  nondimensional  equations 


cosa 

-osina 

osina 

0 

o 

sina 

ocosa 

-ocosa 

0 

d_ 

a 

0 

0 

1 

0 

dx 

0 

-  0 

0 

0 

1 

- 

-  q  - 

28 


-sin0  +  Awsina  +  Atsinat  +  n  -  Acosa 
cos0  -  Awcosa  -  Atcosat  -  Asina 


q 


Vq21* 
L  gr2 


{Ew  +  KAwcosa  -  (l-K)Alcosat) 


CVo  n 

q 

mgrz 


(4.6) 


4.1.2  A  Fictitious  Aircraft 

Level  flight  corresponds  to  y  =  0.  We  assume  that  the  longitudinal  body  reterence  axis 
corresponds  to  the  wing  zero  lift  line  and  that  level  flight  at  nominal  conditions  (Vo.po) 
corresponds  to  a,0  =  ao-  In  this  case,  the  normalized  lift  forces  take  the  form 

Aw  =  fw(a)pu2,  At  =  ft(at)p\)2,  with  fw(0)  =  0,  ft(0)  =  0,  p  =  (p/po)  (4.7) 


The  normalized  drag  force  is  assumed  to  be  of  the  form 
A  =  (a  +  b[fw«x)]2)pu2 

and  the  moment  is  of  the  form 
Iw  =  Gw(ot)p\)2 


(4.8) 


(4.9) 


In  the  following  discussion,  numerical  computations  and  examples  will  be  based  on  the  following 
model  aircraft  characteristics  unless  otherwise  noted. 


p  =  1,  e  =  0,  aw(a)  =  0,  a  =  .1,  b  =.l 


fw  = 


a-2.08(a-ao)3 
ao 


j.  ft  =  e 


/ (a~ao+8e)-3(a-ao+8e 


ao 


^2^  =  300,  =  8 


ao  =  .05,  e'=  .1 


gr 


mgr^ 


(4.10a) 

(4.10b) 

(4.10c) 


4.1.3  Static  and  Dynamic  Stability 


Equation  (4.6)  can  be  written 
J(x)x  =  f(x,u,n) 


(4.11) 
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where  x  denotes  the  state  vector,  u  denotes  the  control  variables  and  |i  denotes  a  designated  set  of 
system  parameters.  A  triple  (x*,u*,|X*)  is  an  equilibrium  point  if 

f(x*,u*,|X*)  =  0  (4.12) 

The  corresponding  perturbation  equations  are 

J(x*)8x  =  ^  (x*,u*,|J*)Sx  +  ^  (x*,u*,p.*)8u  (4. 1 3) 

One  easily  verifies  that  det{  J(x)}=u  for  \)>0,  so  that  under  these  conditions  we  can  rewrite  (4.13) 
in  the  form 

8x  =  A(n*)8x  +  B(fi*)8u  (4. 1 4) 


Example.  Notice  that  k=0,  8e=.0005,n=.l,  11=1,  a=.O495,0=-.O495,q=O  is  an  equilibrium  point,  corresponding  to 
level  flight  at  nominal  velocity.  The  perturbation  equations  are 


8u 

"  -.3960 

-2.949 

-1.0 

0 

d_ 

8a 

-1.980 

-21.80 

0 

1.0 

dt 

80 

0 

0 

0 

1.0 

-  Sq  - 

-  0 

-599.2 

0 

-8.0 

The  dynamical  modes  are 


short  period  X  =  - 14.96 l±j23. 366, 


phugoid  X  =  -  ,13873±j  1.2339, 


-.0053  " 

'  -.0005  " 

-.0300 

I  ' 

.0273 

.0009 

±1 

.0360 

-.8550  - 

-  -.1516  - 

..3723  - 

"  -.2536  " 

-.0081 

±j 

.0043 

-.3847 

-.4069 

.5554  - 

-  -.4182  - 
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In  characterizing  equilibria  for  systems  with  feedback  control  it  is  usual  to  associate  with  (4.1 1)  a 
set  of  outputs  equal  in  number  to  the  control  inputs 

y  =  g(x,u,|i)  (4.15) 

Then  we  obtain  equilibria  by  specifying  p*  and  solving  the  following  equations  for  x*,  u* 


F(x*,u*,p*) 


f(x*,u*,p*) 
g(x*,u*,p*)  . 


(4.16) 


Definition  :  An  equilibrium  point  (x*,u*,p*)  is  regular  if  there  exists  a  neighborhood  of  p* 
on  which  there  exist  unique  functions  x(p),  u(p)  satisfying 


F(x(p),u(p),p)  =0 


with  x*=x(p*),  u*=u(p*). 

Notice  that  the  implicit  function  theorem  implies  that  an  equilibrium  point  is  regular  if 

det[DxFDuF]**0  (4.17) 

Definition  :  An  equilibrium  point  (x*,u*,p*)  is  a  (static)  bifurcation  point  with  respect  to 
F(x,u,p)  if  in  each  neighborhood  of  (x*,u*,p*)  there  exists  (x  1,11141)  and  (x2,U2,p)  with 
(xi,ui)*(x2,U2>  and  F(xi,ui,p)=0,F(x2,U2,p)=0- 

Clearly,  an  equilibrium  point  is  a  bifurcation  point  only  if  it  is  not  regular.  We  can  give  a  useful 
interpretation  to  static  bifurcation  for  systems  defined  by  state  Equation  (4. 1 1)  and  output  Equation 
(4.15).  Let  us  define 

C(p*):=5f  (x*,u*,p*),  D(p*):=g|j  (x*,u*,p*)  (4.18) 

so  that  in  terms  of  perturbation  variables  the  output  equation  becomes 

8y  =  C(p*)5x  +  D(p*)6u  (4.19) 
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With  this  notation  (4.17)  is  equivalent  to 


det 


A(jj*)  B(n*) 
C(^*)  D(^*)  . 


*0 


Thus,  we  have  the  following  conclusion. 


(4.20) 


Theorem  :  An  equilibrium  point  (x*,u*,|i*)  is  a  (static)  bifurcation  point  only  if  the 
linearized  system  (4.14),  (4.19)  has  a  transmission  zero  at  the  origin. 


We  define  stability  of  equilibria  for  parameter  dependent  dynamical  systems  as  follows. 

Definition  :  An  equilibrium  point  (x*,u*,|j.*)  is  (dynamically)  stable  if  it  is  regular  and 
stable  in  the  sense  of  Liapunov. 


Notice  that  this  definition  incorporates  two  essential  elements.  The  equilibrium  point  persists 
under  infinitesimal  variations  of  the  parameter  jj.  (regularity),  and  the  state  trajectories  remain 
bounded  following  sufficiently  small  state  perturbations  (Liapunov  stability). 


It  is  common,  and  useful,  to  distinguish  between  dynamic  and  static  stability.  A  necessary 
condition  for  stability  of  the  equilibrium  point  (x*,u*,p.*)  is  that  the  matrix  A(p*)  has  eigenvalues 
with  nonpositive  real  parts,  so  that  it  must  satisfy 

det{-A(ji*)}  =  (-l)"det{J(x*)}det{|^  (x*,u*,n*)}  2:  0 


where  n=dim(A).  Recall  that  det{J(x)}  =  v  >  0,  so  that  this  condition  reduces  to 
(-l)ndet{^  (x*,u*,p*)}  £  0 


Accordingly,  we  introduce  the  following  notion  of  static  stability. 

Definition  :  An  equilibrium  point  (x*,u*,ii*)  of  (4.11),  (4.15)  is  statically  stable  if  it  is 
regular  and  if  (-lyMet^  (x*,u*,n*)}  >  0. 


Notice  that  static  stability  is  a  necessary  but  not  sufficient  condition  for  stability  Also,  in  the  case 
of  aircraft  longitudinal  dynamics,  as  defined  by  (4.6),  n=4.  Our  definition  of  static  stability  differs 
from  the  conventional  one  (see  Etkin  1 1 1)  only  in  the  fact  that  we  explicitly  include  the  requirement 
that  the  equilibrium  point  be  regular.  It  is  often  the  case  that  aircraft  longitudinal  static  stability 
reduces  approximately  to  the  requirement  that  the  pitch  stiffness  is  negative.  Such  an  example  is 
given  below. 

4.2.  Structured  Model  Uncertainty  of  Aircraft  Longitudinal  Model 

Example  I:  Open  Loop  Properties.  We  give  a  simple  example  which  illustrates  the  importance  of 
center  of  gravity  location  on  aircraft  longitudinal  static  stability.  Consider  the  following  problem. 
With  the  elevator  deflection  angle  8  fixed  and  the  velocity  o  specified,  we  wish  to  determine  values 
of  a,  9  (or,  equivalently,  y)  and  n  which  satisfy  the  equilibrium  equations*,  i  ne  urst  equation 
can  always  be  satisfied  by  choosing 

n  =  sin(a-y)  -  Awsina  -  Atsin(a+S)  +  Acosa  (4. 17) 

Thus,  we  need  only  be  concerned  with  the  determination  of  a  and  0  from  the  remaining  two 
equilibrium  equations 

cos(9)  -  Awcosa  -  Atcos(a+5)  -  Asina  =  0  (4. 1 8a) 

Lw  f  kAwcosoc  -  (l-K)Atcos(a+8)  =  0  (4. 18b) 

Let  us  consider  k  to  be  the  only  adjustable  parameter.  Then  since  a  arid  0  are  the  dependent 
variables  we  have 


DF  =  [DaF  DqF  DkF|  = 


ABO 

.COD. 


A  =  ~  [-Awcoscc-Atcos(a+5)-Asina] 

B  =  -sin9 

C  =  ^  [Lw+KAwcosa-(l-K)Atcos(a+8)] 
D  -  Awcosa  +  Atcos(a+5) 


1  Notice  that  q=0  in  equilibrium. 
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Notice  that  det{Dxf)=-BC=0  only  if  either  0=n7t  for  some  integer  n  or  C=0.  Thus,  a  static 
bifurcation  occurs  only  if  one  of  these  conditions  is  satisfied  simultaneously  with  (4.18).  We  can 
eas;,v  illustrate  the  significance  of  the  case  C=0  as  k  varies.  Equation  (4.18b)  provides  a  relation 
between  the  center  of  gravity  location  (k)  and  the  angle  of  attack  (a).  Figures  4.2  and  4.3  illustrate 
this  relation  for  linear  and  cubic  lift  coefficient  characteristics,  respectively.  Both  curves  reveal 
similar  qualitative  behavior.  There  is  a  critical  eg  location  kc  (and  an  associated  otc,0c) 
corresponding  to  a  local  maximum  and  which  may  be  shown  to  coincide  with  C(ac,0c,Kc)=O.  For 
k>kc  there  are  no  equilibrium  solutions  and  for  k<kc  there  are  two.  In  the  latter  case,  the 
equilibrium  corresponding  to  C<0  (the  one  with  smaller  angle  of  attack)  is  stable  whereas  the  the 
other  equilibrium  corresponds  to  C>0  and  is  unstable. 


eg  Location 


Figure  4.2  Angle  of  attack,  a,  vs.  center  of  gravity  location,  K,  with  d=.4  and  various 
values  of  5  (from  left  to  right  5=-0.01,  0.01,  0.03,  0.04).  Note  that  k  negative  means  that 
the  center  of  gravity  is  forward  of  the  wing  center  of  pressure. 
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eg  Location 


Figure  4.3  Pitch  stiffness,  C,  as  a  function  of  eg  location,  k.  The  same  parameter  values 
as  above. 

Note  that  C  may  be  interpreted  as  the  pitch  stiffness.  It  is  a  commonly  used  indicator  of  aircraft 
longitudinal  static  stability  [1-4].  In  the  aviation  community  static  stability  is  distinct  from  dynamic 
stability.  An  equilibrium  point  is  said  to  be  statically  stable  if  C>0,  statically  unstable  if  C<0  and  to 
have  neutral  static  stability  if  C=0.  In  the  proceeding  example  neutral  static  stability  corresponds  to 
a  parameter  value  at  which  the  equilibrium  point  is  not  regular  -  indeed,  it  corresponds  to  a 
bifurcation  point. 

It  should  be  emphasized  that  this  static  instability  is  distinct  from  aerodynamic  stall  phenomena. 
The  cubic  lift  coefficients  used  in  generating  Figure  4.3  are  illustrated  in  Figure  4.4.  It  is  readily 
observed  that  the  critical  angle  of  attack  in  Figure  4.3  is  well  below  the  stall  angle  of  attack  for 
either  the  wing  or  the  tail. 
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Figure  4.4.  Normalized  lift  coefficients  fw  (solid)  and  fye'  (dashed)  plotted  as  functions  of 
a. 

In  the  following  Figures  4.6  through  4.1 1  we  illustrate  how  velocity  affects  the  equilibrium  values 
of  pitch  angle  and  angle  of  attack. 


Figure  4.5a  Pitch  attitude, 0,  vs.  eg  location,  K,  with  u=.4,  and  from  left  to  right  S=-0.01, 
0.01,  0.03,  0.04,  0.05,  0.06. 


36 


Pitch  Attitude 


Figure  4.5b  Angle  of  attack, c 
0.01,  0.01,  0.03,  0.04,  0.05,  I 


eg 

Figure  4.6a  Pitch  attitude,©, 
0.01.  0.01.  0.03.  0.04.  0.05.  I 


Pitch  Attitude  Angle  of  Attack 


eg  Location 


Figure  4.6b  Angle  of  attack, a,  vs.  eg  location,  k,  with  1o=.41,  and  from  left  to  right  6 
0.01,  0.01,  0.03,  0.04,  0.05,  0.06. 


Figure  4.7a  Pitch  attitude,©,  vs.  eg  location,  k,  with  u=.42,  and  from  left  to  right  8 


0.01,  0.01,  0.03,  0.04,  0.05,  0.06. 
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Attitude 


eg  Location 


Figure  4.7b  Angle  of  attack.oc,  vs.  eg  location,  K,  with  0=42,  and  from  left  to  right  5 
0.01,  0.01,  0.03,  0.04,  0.05,  0.06. 


Figure  4.8a  Pitch  attitude,©,  vs.  eg  location,  K,  with  v=.45,  and  from  left  to  right  5 
0  01,  0.01,  0.03,  0.04,  0.05,  0.06. 


Attitude  Angle  of  Attack 


Pitch  Attitude 


eg  Location 


Figure  4.9b  Angle  of  attack.a,  vs.  eg  location,  k,  with  \)=.6,  and  from  left  to  right  6 
0.01,  0.01,  0.03,  0.04,  0.05,  0.06. 


Figure  4.10a  Pitch  attitude,9,  vs.  eg  location,  K,  with  o=.8,  and  from  left  to  right  o 
0.01,  0.01,  0.03,  0.04,  0.05,  0.06. 


Pitch  Attitude 


I 

I 


Figure  4.11a  Pitch  attitude, 0,  vs.  eg  location,  K,  with  v=l.,  and  from  left  to  right  8 
0.01,  0.01,  0.03,  0.04,  0.05,  0.06. 
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Angle  of  Attack 


Angle  of  Attack  Pitch 
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eg  Location 

Figure  4.12b  Pitch  Stiffness,  C,  vs.  eg  location,  K,  with  o=.42,  and  elevator  deflection 
angle  8=  0.03. 
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Figure  4.12c  Angle  of  attack, a,  vs.  eg  location,  k,  with  u=.42,  and  elevator  deflection 
angle  8=  0.03. 
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Figure  4. 13a  Pitch  attitude,0,  vs.  eg  location,  k,  with  v>=l and  elevator  deflection  angle 
6=  0.03. 
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Figure  4.13b  Pitch  Stiffness,  C,  vs.  eg  location,  tc,  with  \)=1.,  and  elevator  deflection 
angle  5=  0.03. 


Figure  4.13c  Angle  of  attack, a,  vs.  eg  location,  k,  with  t>=l.,  and  elevator  deflection 
angle  8=  0.03. 

Example  II:  Closed  Loop  Properties.  A  somewhat  more  pertinent  example  with  respect  to  control 
system  design  is  the  following.  Once  again  consider  the  longitudinal  dynamics  defined  by 
equation  (4.6).  It  is  desired  to  regulate  the  velocity  and  flight  path  angle  y  by  adjusting  the 
elevator  deflection  angle  and  thrust  8,  FI.  Thus,  we  define  the  output  equations 


'  yi  ' 

["  v-v*  " 

v-x>* 

-  y2  - 

L  y-y*  - 

.  a-Q-y*  . 

Given  the  desired  flight  path  parameters  v,  y  we  wish  to  determine  values  of  a,  8,  n  which  satisfy 
the  equilibrium  equations  (4.16),  with 


f(x,u,(l) 


-sin0  +  Awsina  +  Atsinat  +  II  -  Acosa 
cos8  -  Awcosa  -  Atcosat  -  Asina 

q 

V02l*  CVn 

{Lw  +  KAwcosa  -  (l-K)Atcosat)  -  — H,  q 
gX  mgrz 


(4.19a) 


g(x,u,n) 


d-d* 
a-0-y*  . 


(4.19b) 
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Once  again  n  is  directly  determined  from  (4. 17)  and  we  need  only  be  concerned  with  the  solution 
properties  of  the  pair  of  equations 


cos(a-y)  -  Awcosa  -  AtCos(a+5)  -  Asina  =  0  (4.20a) 

Lw  +  KAwcosa  -  (l-K)Atcos(a+5)  =  0  (4.20b) 


Figures  4.14  a&b  illustrate  the  equilibrium  values  of  elevator  deflection,  5,  and  angle  of  attack,  a, 
as  a  function  of  eg  location,  k,  at  cruise  conditions  u=l,  Y=0. 


Figure  4. 14a  Elevator  deflection,  5,  as  a  function  of  eg  location,  k,  witht>=l ,  y=0. 


eg  Location 

Figure  4. 14b  Angle  of  attack,  a,  as  a  function  of  eg  location,  tc,  with\)=l ,  y=0. 
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In  Figures  4.15a&b  we  illustrate  the  eigenvalue  locations  as  a  function  of  eg  location,  (Kk<.4. 
Figure  4.15a  clearly  shows  the  expected  short  period  and  phugoid  branches.  A  loss  of  static 
stability  occurs  at  approximately  k=.12.  Interestingly  enough,  although  the  Jacobian 
Dxf(x*,u*,|i*)  is  singular,  this  does  not  correspond  to  a  bifurcation  point  because  the  Jacobian 
Dx,pF(x*,u*,|i*)  is  not. 


Figure  4.15a  Open  loop  eigenvalue  locations  as  a  function  of  eg  location,  k,  at  cruise 
conditions. 


Figure  4.15b  Near  origin  blowup  of  open  loop  eigenvalue  locations  as  a  function  of  eg 
location,  K,  at  cruise  conditions.  This  figure  clearly  illustrates  occurrence  of  a  divergence 
instability  as  the  eg  moves  toward  the  rear  of  the  aircraft  and  one  of  the  phugoid  roots 
crosses  into  the  right  half  plane.  This  occ  "s  at  approximately  tc=.12. 
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4.3.  Worst  Case  Design  Analysis  for  Aircraft  with  Relaxed  Static  Stability 
The  linear  dynamic  model  for  the  aircraft  described  in  the  previous  section  taken  about  the  nominal 
operating  condition  k=0  is  statically  and  dynamically  stable.  The  linear  perturbation  equations  can 
be  written  in  the  state  space  form  as 
x  =  Akx  +  Bku 
y  =  cKx 

with  state  x=[t),y,0,q]t,  control  u=[n,8e]t.  For  K=0  we  obtain  the  equilibrium  conditions  for 
level  flight  in  terms  of  the  state  vector. 


1. 

0.0497523 

0.0497523 

0 


and  the  control 
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eq 


-  0.0002477 ' 
.  0.0996290 . 


The  linear  perturbation  model  coefficients  are  then  obtained: 
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The  model  exhibits  the  characteristic  pair  of  underdamped  short  period  and  phugoid  dynamic 
modes; 

Open  loop  poles  at  k=0 


-  11.587041  +22.020408i 
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-  11.587041  -22.020408i 

-  0.0622142  +  1.386446H 

-  0.0622142  -  1.3864461i 


The  choice  of  outputs  defined  by 


C  = 


'  1.  0.  0.  o.- 
.0.  -1.  1.  0.. 


gives  a  real  transmission  zeros  at  77.196225  and  -  78.196225.  The  right  half  plane  transmission 
zero  which  will  ultimately  limit  the  bandwidth  of  the  closed  loop  system  for  robust  stabilization. 
The  pole  zero  plot  is  shown  in  Figure  4.16. 


X  poles 

•  zeros  basile 


Figure  4. 16:  Pole-zero  plot  for  nominal  statically  stable  aircraft. 

If  the  CG  shifts  aft  by  .12  relative  to  the  effective  nominal  moment  arm  for  tail  controlled  pitch 
motion  the  aircraft  becomes  both  statically  and  dynamically  unstable.  The  resulting  pole- zero  plot 
is  shown  in  Figure  4.17,  although  the  equilibria  shift  only  slightly  to  : 
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Figure  4.17:  Pole-zero  plot  for  statically  unstable  aircraft  configuration. 

To  procede  with  the  minimax  (worst  case)  design  we  first  determine  the  minimum  achievable 
stability  margin  based  on  the  optimization  procedure  and  its  dependence  on  the  unknown  parameter 
0  <k<0.12.  This  dependence  is  shown  in  Figure  4.18.  The  optimal  geometric  stability  margin  is 
dependent  on  the  CG  shift  as  shown.  For  this  range  of  the  independent  parameter  the  worst  case 
condition  is  the  aft  most  location  of  the  CG  as  expected. 
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Figure  4.18:  Optimal  Stability  Margin  Dependence  on  k. 

The  optimal  robust  geometric  stability  margin  can  be  given  the  interpretation  of  providing  a  bound 
on  the  allowable  perturbation  to  a  normalized  left  coprime  factorization  of  the  plant.  One  can 
therefore  ask  if  the  unstructured  uncertainty  condition  for  robust  synthesis  can  provide  conditions 
for  stabilization  of  the  flight  dynamics  model  with  aftward  shift  in  the  CG  of  0.12.  From  the 
above  the  optimal  stability  margin  variation  is  continuous  and  monotonically  decreasing  with  the  aft 
shift.  Taking  the  aft-most  location  as  worst  case  we  can  compute  a  maximum  additive  plant 
perturbation  directly  in  terms  of  the  normalized  left  coprime  factors  of  the  plant  for  the  nominal 
case 

G0(s)  =  Mq(s)N0(s) 


and  for  the  case  at  k=0. 12; 

Gk(s)  =  s) 

then  we  can  obtain  the  required  bounds  by  plotting  the  singular  values  of 

[  M0  NQ]  (jco)  -  [  5^  (jto) 

as  shown  in  Figure  4.19. 
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Figure  4.19:  Singular  Values  of  NCP  Factor  Perturbation 

Thus  using  the  unstructured  model  uncertainty  description  obtained  from  the  geometric  stability 
margin  and  the  expected  perturbation  to  the  NLCP  factors  we  see  that  the  required  bound  is 

esm=0-7619.  This  means  that  to  design  a  stabilizing  compensator  using  a  totally  unstructured 

requirement  is  so  conservative  in  this  case  that  it  cannot  in  fact  be  obtained!  Indeed,  from  the 
above  plot  of  optimal  stability  margin  vs.  CG  shift  we  see  that  no  operating  point  obtains  an 
optimal  geometric  stability  margin  in  excess  of  0.56.  This  is  typical  of  unstructured  uncertainty 
modeling  and  robust  control  synthesis.  However,  it  is  significant  to  note  that  using  the  geometric 
stability  margin  and  the  optimal  robust  stabilization  problem  a  unique,  quantitative  criterion  for 
robust  stabilization  is  obtained  which  addresses  the  problem  without  examining  frequency 
dependent  data. 

To  illustrate  the  options  for  design  using  the  optimal  robust  stabilization  formulation  we  first 
consider  scaling  and  its  effect  on  the  closed  loop  design.  We  first  consider  the  design  for  a  fixed, 
nominal  operating  point. 

The  MIMO  open  loop  system  has  gains  displayed  as  frequency  dependent  singular  value  plot  as 
shown  in  Figure  4.20. 
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Figure  4.20:  Singular  value  Loop  Gains  for  nominal  Plant  Transfer  Function. 


Given  the  ultimate  bandwidth  limitation  imposed  by  the  right  half  plane  zero  at  -77  we  choose  a 
nominal  crossover  frequency  of  <0=10.  Rescaling  the  inputs  and  outputs  proceeds  first  by 
consideration  of  the  relative  importance  of  changes  in  each  channel.  For  the  nondimensional  model 
the  inputs  consist  of  thrust  and  deflection  of  the  tail  control  surface  and  the  outputs  include  velocity 
and  flight  path  angle.  Scaling  to  obtain  equalization  of  each  channel  follows  from  the  choice  of 
nominal  value  and  relevant  changes  as: 

nominal  relative 


r°l  r  0.1  n 

L  0  J  L  0.5236  J 

r 1 ]  r  01  ] 

L  0  J  L  .0873  J 


Let  the  input/output  scaling  be  given  by  the  diagonal  matrices  as : 


r  .1  0  1 

U  =  SiI=Lo  .5236  J® 

r  .05  On 
y*S°H  0  .0873  j 


then  the  transfer  function  becomes  G(s)  =  S'J G(s)Sf . 
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Finally,  to  specify  desired  crossover  frequency  we  reduce  the  effective  plant  gain  in  each  channel 
so  that  the  maximum  singular  value  at  desired  crossover  satisfies  .  The  resulting  scaled  plant  has 
gains  as  shown  in  Figure  4.21  for  the  worst  case  parameter  value  of  k=.12. 


Figure  4.21:  Singular  Value  Loop  Gains  of  Scaled  Plant 


The  optimal  geometric  stability  margin  for  the  scaled  plant  is  now  computed  as  esm=.2772.  We 

note  that  the  procedure  described  previously  suggests  that  the  system  has  compatible  loop  shape. 
We  can  now  obtain  a  suboptimal  controller  for  any  desired  level  of  stability  robustness  as 
y>l/.56i=1.783.  We  take  the  suboptimal  level  y=1.8  and  obtain  by  the  (2.19)-(2.22)  the 


controller  realization  as: 

- -137.81869  -159.443  157.59503  0.0071676  ” 

A  = 

202.57323  240.27344  -262.55426  0.9675697 

c 

687.31722  884.00473  -884.00473  1 

-  4043.8865  4666.5232  -5317.9173  -10.610107- 

B  = 

”  -6.694682  13.840263  ” 

10.200762  -22.84293 

c 

34.365861  -77.173613 

-  196.7627  -  440.94193- 

f  4.7375856  9.371697  -11.645754  -0.41910801 

C  = 

c  L  18.620063  -0.4424553  -0.2898523  -0.0357794. 


To  illustrate  that  the  resulting  suboptimal  controller  achieves  the  required  stability  margin  over  a 
broad  range  of  frequencies  we  plot  the  singular  values  of 

O(jto)  =  [MQ  -  RPl(j(0) 

where  K(s)  =  P(s)Q_1(s)  is  a  normalized  right  coprime  factorization  for  the  achieved  controller 
K(s)  =  QlsI-AJ^Bc  (see  Figure  4.22).  To  highlight  the  broadband  solution  obtained  from  this 
problem  we  plot  separately  the  maximum  singular  value  of  F  (i.e.  the  geometric  stability  margin)  in 
Fig.  4.22a  and  minimum  singular  value  in  Fig.  4.22b. 
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Figure  4.22a:  Maximum  Singular  Value  of  <1> 
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Figure  4.22b:  Minimum  Singular  Value  of  <1>. 

Finally,  the  achieved  loop  gains  are  displayed  by  plotting  the  singular  values  of  K(s)G(s)  in  Figure 
4.23. 


Figure  4.23:  Achieved  loop  transmission  gains. 


For  direct  comparison  we  superimpose  the  desired  loop  shape  with  that  achieved  from  the 
suboptimai  compensator  for  the  scaled  plant  response. 
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Figure  4.24:  Achieved/desired  loop  transmission  gains 


4.4.  Performance  Considerations  for  Aircraft  Flight  Control 

In  this  problem  we  introduce  two  modifications  of  the  previous  problem.  First  to  illustrate  design 
for  performance  robustness  we  introduce  a  requirement  for  integral  action  in  each  channel. 
Second,  to  test  the  worst  case  parametric  design  problem  we  perform  the  design  for  the  worst  case 
operating  point  at  k=.12.  We  also  assume  there  is  a  dominant  source  of  modeling  error  which  is 
characterized  as  relative  error  in  terms  of  bounds  on  multiplicative  perturbation  at  the  system 
output.  We  assume  this  is  given  by  a  specification  that  the  model  is  to  be  considerd  accurate  to 
10%  for  frequencies  less  than  2.  Beyond  co>2  the  model  relative  error  increases  at  20dB/decade. 
Thus  the  relative  bound  lm>HDmll  is  greater  than  1  for  to>50. 

To  motivate  the  choice  of  weighting  functions  we  consider  the  requirement  for  steady  error 
reduction  and  consequences  for  integral  action  in  each  loop.  From  the  discussion  in  Blight, 
Gangsaas  and  Richardson  [1]  we  see  such  requirements  are  apparent  in  LQG/LTR  type  designs  but 
the  application  is  much  more  difficult  to  implement  than  using  the  present  method.  To  balance  the 
requirements  for  integral  action  with  loop  crossover  at  co=10  we  performed  a  simple  tradeoff  with 
PI  type  weighting  in  each  channel  to  obtain  the  desired  loop  shape.  The  weighting  considered  has 
the  form: 


Wi(s)  = 


— +1 
tjS 

0 


0 


TjS 
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A  value  of  \  =.l  was  chosen. 


It  became  apparent  that  from  the  perspective  of  multiplicative  model  errors  at  the  outputs  the  high 
frequency  attenuation  needed  to  be  increased.  We  therefore  added  first  order  rolloff  in  terms  of  the 
leadylag  weights: 


W2(s)  = 


TzS+1 

TpS+1 

0 


0 

TzS+1 

TpS+1 


where  values  of  8  and  80.  were  chosen  for  the  lead  zero,  resp.  lag  pole.  This  permitted  increased 
rolloff  in  the  region  after  cutoff  without  serious  attenuation  near  the  desired  cutoff  frequency.  The 
shaped  loop  gain  in  displayed  in  Figure  4.25. 


frequency 

Figure  4.25:  Desired  loop  shape  with  weightings. 


The  shaped  transfer  function  dynamics  are  summarized  in  terms  of  the  transmission  poles  and 
zeros  in  Table  4.1. 


Table  4.1:  Shaped  Transmission  Poles  and  Zeros 


poles 

zeros 

-  80. 

77.196225 

-  80. 

-  78.196225 

0 

-  8. 

0 

-  8. 

29.556023 

-  0.7 

S.79 12262 

-0.7 

-  0.2004232  +  1.4435674i 
-0.2004232-  1.4435674i 

For  this  loop  shape  we  compute  the  optimal  stability  margin  as  £sm=0.2772.  Then  we  choose  the 
desired  suboptimal  level  for  controller  synthesis  as  y>l/.2772.  Here  we  take  y=3.65  and 
synthesize  the  suboptimal  robust  control  for  the  weighted  plant.  Incorporating  the  weighting 
functions  into  the  controller  K(s)  we  obtain  the  achieved  loop  gains  as  displayed  in  Figure  4.26.  It 
can  be  seen  that  the  specifications  given  in  terms  of  relative  model  uncertainty  at  the  outputs 
(multiplicative  errors)  together  with  performance  specification  for  integral  action  and  cross  over 
frequency  at  least  (D=10  have  been  achieved. 


fjpqwry 

Figure  4.26:  Acheived  Loop  Gains  for  shaped  worst  case  design. 


To  check  the  assumptions  implicit  in  the  minimax/maximin  equivalence  for  parametric  robustness 
we  evaluate  the  stability  margin  achieved  by  the  worst  case  design  at  the  nominal,  stable 
configuration  of  k=0  (the  other  extreme  of  the  parameter  variation  in  this  case).  The  resulting 
stability  margin  achieved  is  esm=0.2583  as  seen  from  Figure  4.27. 
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Figure  4.27:  Achieved  Stability  Margin  @k=0. 

Finally,  to  illustrate  the  effect  of  the  worst  case  design  on  the  nominal,  stabel  aircraft  configuration 
we  display  the  achieved  loop  gains  for  this  case  in  Figure  4.28. 
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Figure  4.28:  Achieved  Loop  Gains  for  stable  configuration. 


The  closed  loop  system  poles  for  the  nominal  stable  configuration  with  the  control  as  designed 
above  can  be  computed  as: 


Closed  Loop  Poles  for  nominal  (stable)  aircraft 

-  1819.4462 

-  270.14803  +  29.621 8 17i 

-  270.14803  -  29.621817i 

-  78.048361 

-  7.222071  +  25.3900961 
7.222071  -  25.390096i 

-  12.391449 

-  7.4595557 

-  2.0803007 

-  2.0285601 

-  0.7878369 

-  0.6121988 

5.  Conclusions  and  Directions  for  Future  Research 

5.1.  Nonlinear  Control  Design  for  Robust  Flight  Control 

In  order  to  achieve  enhanced  maneuverability  and  efficiency,  future  aircraft  will  operate  close  to  or 
even  beyond  open  loop  stability  boundaries.  For  example,  reduction  of  horizontal  tail  size  in  order 
to  achieve  reduced  fuel  consumption  results  in  loss  of  longitudinal  static  stability  for  sufficiently  aft 
c.g.  locations  [  1 ).  Fighter  aircraft  may  operate  at  high  angle  of  attack  or  at  high  roll  rate  where 
nonlinear  effects  cause  loss  of  stability  16].  Such  aircraft  require  augmentation  by  automatic  flight 
control  systems  which  induce  the  desired  handling  qualitites  over  the  full  range  of  flight 
conditions. 

Feedback  controllers  for  these  applications  are  typically  conceived  to  be  linear,  perhaps  with  some 
form  of  gain  scheduling.  However,  when  operating  near  stability  boundaries  the  system  dynamics 
are  nonlinear  in  an  essential  way.  It  is  well  known,  for  instance,  that  a  linear  perturbation  model  is 
a  reliable  indicator  of  stability  only  if  none  of  its  eigenvalues  lie  on  the  imaginary  axis  (9). 
Otherwise  nonlinear  effects  are  crucial.  Some  recent  studies  in  flight  mechanics  characterize 
aircraft  loss  of  stability  in  terms  of  elementary  local  bifurcations  1 10-11],  Thus,  divergence 
instability  is  typically  associated  with  a  saddle-node  b:rurcation  and  a  change  in  the  equilibrium 
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point  structure  while  flutter  is  associated  with  a  Hopf  bifurcation  and  the  appearance  of  a  limit 
cycle.  Such  phenomenon,  of  course,  are  fundamentally  nonlinear. 

Unfortunately,  these  studies  deal  almost  exclusively  with  open  loop  dynamics  under  parameter 
variation  and  very  little  general  theory  is  availiable  concerning  the  design  of  feedback  controls  near 
bifurcation  points.  We  can,  however,  draw  some  obvious  inferences  although  our  remarks  will  be 
confined  to  static  instability  (divergence).  It  is  known  that  linear  feedback  can  be  used  to  stabilize  a 
divergence  instability.  McRuer  et  al  |4|  note  that  in  aircraft  such  stabilization  requires  high  gain, 
wide  bandwidth  control.  Such  controllers  may  be  sensitive  to  actuator  saturation  and  excitation  of 
high  frequency  parasitics.  It  is  our  view  that  these  problems  assume  exagerated  proportions 
because  of  the  attempt  to  force  a  linear  solution  on  an  intrinsically  nonlinear  problem. 

Note  also  that  when  the  system  operates  near  a  saddle-node  bifurcation  point,  the  stable 
equilibrium  point  is  necessarily  close  to  the  boundary  of  the  domain  of  attraction  (the  neighboring 
unstable  equilibrium  is  on  the  boundary.)  It  follows  that  the  system  may  be  unacceptably  sensitive 
to  external  disturbances.  Performance  can  sometimes  be  dramatically  improved  by  nonlinear 
feedback.  In  fact,  the  substantial  improvement  obtained  by  Garrard  and  Jordan  |6|  in  recovery 
from  stall  by  using  a  nonlinear  feedback  can  be  explained  by  this  observation. 

Techniques  for  the  design  of  nonlinear  feedback  control  systems  are  still  very  limited  and  tend  to 
be  tailored  to  specific  situations.  The  most  promising  approaches  appear  to  be  those  associated 
with  methods  of  exact  linearization.  This  procedure  is  based  on  some  early  work  of  Krener  1 44 1 
and  others  demonstrating  that  a  large  class  of  state  equation  models  for  nonlinear  dynamical 
systems  can  be  exactly  linearized  by  a  combination  of  nonlinear  stale  feedback  and  a  nonlinear 
transformation  of  state  coordinates.  Once  the  system  is  linearized  all  methods  of  control  system 
design  for  linear  dynamics  become  applicable.  The  inverse  transformation  is  then  implemented  in 
conjunction  with  the  linear  compensator  so  that  the  overall  controller  is  nonlinear.  Meyer  and  his 
coworkers  have  articulated  the  application  of  these  ideas  to  certain  problems  in  flight  control  as  an 
alternative  to  gain  scheduling. 

Another  approach  termed  input output  linearization  has  evolved  from  work  of  Hirschorn  |4.5|.  The 
notion  of  input-output  linearization  as  developed  by  Kravaris  and  Chung  |46| )  is  based  on  the  so- 
ealled  Myrtles  Isidori  canonical  form  1 47 1  for  nonlinear  systems.  This  approach  produces  a  linear 
input  output  model  and  may  have  certain  advantages  especially  when  frequency  domain  design 
methods  are  to  be  used  for  compensator  design.  The  attractiveness  of  such  methods  is  further 
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enhanced  by  recent  frequency  domain  formulations  of  Hopf  bifurcation  analysis  [12]  with 
application  to  multivariable  nonlinear  control  system  analysis  [54], 


Although  these  methods  may  be  formally  applied  at  or  near  bifurcation  points  such  applications 
have  not  been  studied.  It  is  interesting  to  note  that  another  approach  to  linearization  which  has 
generated  some  interest  in  recent  years,  the  extended  linearization  method  of  Rugh  [59]  fails  in  an 
obvious  way  at  bifurcation  points. 

The  essentials  of  the  approach  are  most  easily  understood  in  terms  of  the  single  input  single  output 
problem.  Consider  a  nonlinear  dynamical  system  in  the  form 


X  =  f(x)  +  g(x)u 

(5.5a) 

y  =  h(x) 

(5.5b) 

where  f,  g  are  smooth  vector  fields  on  Rn  and  h  is  a  smooth  function  mapping  Rn  — »  R.  Now, 
differentiate  (5.5b)  to  obtain 

y  =  ^(f(x)  +  g(x)u)  (5.6) 

If  the  scalar  coefficient  of  u  is  zero,  we  differentiate  (5.6)  and  continue  in  this  way  until  a  nonzero 
coefficient  first  appears. 

This  process  can  be  succinctly  described  by  introducing  some  conventional  notation  of  differential 
geometry  and  analytical  mechanics.  We  need  only  the  concept  of  Lie  derivative.  First,  the  Lie 
(directional)  derivative  of  the  scalar  function  h  with  respect  to  the  vector  field  f  is  defined  as 

Lrfh)  =  <dh,f>  =  gj|f(x)  (5.7) 

Since  the  Lie  derivative  is  itself  a  scalar  function  on  Rn,  higher  order  derivatives  may  be 
successively  defined 

L* (h)  =  LKL^(h))  =  <dL^(h),  f )  (5.8) 

Now,  (5.6)  can  be  written 
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y  =  <dh,  f  )  +  (dh,  g  )u  =  Lf(h)  +  Lg(h)u  (5.9) 

If  Lg(h)  =  0,  then  differentiate  (5.9)  to  obtain 

y  =  (dLf(h),f )  +  (dLf(h),g  >u  =  L*(h)  +  Lg(Lf<h))u  (5. 10) 

If  Lg(Lf(h))  =0,  then  differentiate  (5.10).  If  Lg(L^_1(h))  =  0  for  k  =  1,..,  r-1,  but  Lg(Lrf  1  (h))  *  0, 
then  the  process  ends  with 

^  =  Lrf(h)  +  Lg(L^,(h))u 

The  number  r  is  called  the  “relative  degree”  or  “characteristic  number”  of  (5.5) 
define  the  coordinates  zg  Rr 
zk  =  Lj'1(h),  k  =  l,..,r 

then  Equations  (5.9)  through  (5.1 1)  can  be  written 
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_  a(x)  +  p(x)u  _ 

where 

a(x)  =  L^(h),  and  p(x)  =  Lg(L^.'1(h)) 

Formally,  the  system  (5.12)  may  be  linearized  and  simultaneously  stabilized  in  the  input-output 
sense  as  follows.  Choose  the  control 

u  =  <D(x,v)  =  (v  -  a(x))/p(x)  (5.13) 


(5.11) 

.  Note  that  if  we 


where 


r-1 


O(x)  =XPkL?(h)  +  Lj(h), 

k=0 


p(x)  =  Lg(L[’1(h)) 
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and  (3k,  k  =  0,..,  r-1  are  real  numbers  corresponding  to  the  coefficients  of  any  Hurwitz 
polynomial.  Then  the  relationship  between  the  new  control  input  v  and  the  output  y  is  defined  by 
the  linear  r-dimensional  (completely  controllable  and  observable,  companion  form)  system 


(5.14a) 


y  =  [  1  0  .  .  0]z  (5.14b) 

These  calculations  are  readily  extended  to  MIMO  case  [55,  56]. 

The  importance  of  this  construction  is  that  linear  dynamic  compensators  may  now  be  designed  for 
the  resulting  linear  system  using  any  applicable  linear  method.  The  exact  linearization  is  global  if 
the  transformations  admit  unbounded  control  u;  otherwise,  the  linearization  is  local.  However,  it 
should  be  emphasized  that  unlike  traditional  linearization  about  a  fixed  operating  point  or  trim 
condition  (based  formally  on  Taylor  series  expansion)  the  approach  just  described  linearizes  the 
system  dynamics  about  a  nominal  model.  Open  questions  in  the  available  theory  for  such 
linearization  includes  the  effect  of  such  control  constraints  and  parametric  uncertainty. 

5.2.  Linear  Mode  Design  and  Critical  Nonlinear  Dynamics 

For  the  limited  scope  of  the  phase  1  study  we  have  focused  on  computational  studies  for  relaxed 
static  stability  aircraft  where  the  parametric  uncertainty  is  related  to  the  shift  of  CG  during  flight 
operations.  The  underlying  nonlinear  dynamics  illustrate  the  role  of  parameters  in  determining  the 
dynamic  properties  of  the  resulting  family  of  linear  perturbation  models.  What  is  evident  from  our 
analysis  and  numerical  computations  so  far  is  the  significance  of  the  transmission  zeros  to  this 
parameter.  Certainly  in  the  region  which  we  considered  for  variation  of  K  the  transmission  zero  are 
relatively  insensitive  by  comparison  with  the  system  poles.  However,  this  does  change  as  the 
variation  of  k  is  extended  to  approach  the  static  bifurcation  point  for  the  system  equilibria.  It  has 
been  shown  that  in  feedback  systems  static  bifurcation  is  associated  with  the  passing  of  a  real  zero 
through  the  origin.  This  differs  from  the  conventional  case  in  dynamical  sytem  theory  where  static 
bifurcation  is  associated  with  a  real  eigenvalue  passing  through  the  origin.  In  this  region  we  expect 
the  sensitivity  of  the  system  optimal  stability  margin  to  be  significant,  indicating  a  difficult  control 
design  problem.  It  is  our  conjecture  that  part  of  the  limitation  evidenced  in  such  regions  comes 
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from  the  linear  mcxle  design  approach  where  we  are  essentially  trying  to  design  a  single  fixed  (and 
therefore  robust)  linear  compensator  for  the  linear  plant.  The  requirement  for  stability  of  the  linear 
mode  design  has  the  implicit  objective  to  maintain  the  system  state  near  the  desired  equilibrium  or 
trim  condition.  Near  a  bifurcation  of  the  equilibria  under  parametric  variation  we  expect  the 
domain  of  attraction  to  change  depending  on  the  type  of  bifurcation.  Clearly,  the  underlying 
nonlinear  dynamics  will  be  significant  for  control  of  the  system  dynamics  in  such  operating 
regimes. 

Such  considerations  have  lead  engineers  to  develop  various  methods  for  gain  scheduling-most  of 
which  are  ad  hoc.  On  the  other  hand  the  previous  discussion  on  feedback  linearizing  control  for 
certain  nonlinear  systems  provides  a  complete  analytical  basis  for  control  design  based  on 
linearization  about  a  nominal  system  model  rather  than  about  an  equilibrium  operating  point  in  the 
system  state  space.  Despite  the  abstract  basis  for  the  theory  of  feedback  linearization,  its  practical 
application  provides  an  analytical  basis  for  gain  scheduling  the  linear  mode  design.  Both  methods 
are  essentially  model-based  control  design  schemes  and  therefore  robustness  considerations  are 
important  in  applications.  Flight  control  systems  often  employ  gain  scheduling  with  respect  to 
various  measurable  system  parameters  (e.g.  velocities,  angle  of  attack,  etc.)  for  operation  in 
nonlinear  regimes.  It  is  our  view  that  methods  for  integrated  design  of  multiloop  controllers  based 
on  methods  related  to  feedback  linearization  may  provide  enhanced  capabilities  for  such  high 
performance  flight  control  problems  as  recovery  from  stall.  In  the  Phase  2  effort  we  will  propose 
to  study  methods  for  robust  design  and  implementation  of  such  nonlinear  control  laws. 

5.3.  Research  Directions  for  multiparameter  worst  case  design 

The  worst  case  design  procedure  based  on  the  maximin  design  is  computationally  feasible  if  the 
number  of  real  parameters  is  small  since  the  evaluation  of  the  optimal  controller  dependence  on  a  is 
simple.  It  is  therefore  suggested  that  in  applications  the  plant  be  modeled  with  the  less  significant 
real  parametric  uncertainty  embedded  in  the  unstructered  uncertainty  bounds.  This  is  the  practical 
aspect  of  the  methods  we  will  propose  to  investigate  and  develop  in  Phase  2:  that  we  seek  a 
pragmatic  approach  to  implement  the  available  methods  in  optimization  based  approachs  to  robust 
control  synthesis  for  structured  uncertainty  together  with  more  well  developed  methods  for 
structured  uncertainty.  We  believe  that  this  blend  can  be  developed  with  respect  to  a  significant 
class  of  flight  control  problems  where  dynamic  degrees  of  freedom  are  relativly  small. 

The  focus  of  the  Phase  1  study  has  been  on  developing  a  method  for  robust  stabilization  of  models 
subject  to  combined  real  parametric  uncertainty  and  parasitic  (unmodeled)  dynamics.  Our  efforts 
were  directed  toward  eliminating  the  source  of  conservatism  abundant  in  such  design  methods 
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while  retaining  the  essential  feature  of  robust  design  in  the  engineering  context.  We  feel  that  the 
optimal  worst  case  design  methods  described  above  addresses  these  concerns  in  principle.  The 
result  however,  is  a  computationally  difficult  problem  in  the  general  case.  Our  approach  in  the 
benchmark  designs  has  been  pragmatic.  We  also  believe  that  the  essential  features  of  the  nonlinear 
dynamics  of  flight  control  designs  will  serve  to  restrict  consideration  of  the  variation  of  physically 
based  real  parameters  which  effect  the  aircraft  dynamics.  Nevertheless,  we  feel  that  a  focus  of  the 
Phase  2  proposal  should  be  in  the  area  of  algorithm  development  for  the  minimax  (worst  case) 
design  for  realistic  combination  of  structured  uncertainty  (arising  from  physical  parameters)  and 
unstructured  uncertainty  (arising  from  unmodeled  dynamics). 

We  also  feel  that  the  entropy  interpretation  of  the  results  of  Bernstein  and  Haddad  —and  therefore 
the  results  obtained  in  Yeh  et  al  [41]  —  may  offer  additional  perspectives  on  how  inherent 
consc-vatism  of  such  designs  can  be  reduced.  The  analysis  achieved  to  date  represents  only  a 
special  case  for  these  problems. 

The  loop  shaping  procedure  has  been  shown  to  be  effective  with  many  desireable  properties. 
However,  the  choice  of  loop  shape  has  not  been  adequately  addressed.  In  particular,  the 
modification  of  loop  shape  to  incorporate  knowledge  about  system  uncertainties,  the  question  of 
how  loop  shape  affects  Emax>  and  the  role  and  use  of  loop  shaping  for  decoupling  and  scaling  are 
all  potentially  fruitful  directions  of  inquiry. 
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ABSTRACT  -  A  two  stage  Hx  based  design  procedure  has 
been  described  which  uses  a  normalized  coprime  factor  approach 
to  robust  stabilization  of  linear  systems.  A  loop-shaping  pro¬ 
cedure  is  also  incorporated  to  allow  the  specification  of  perfor¬ 
mance  characteristics.  Theoretical  justification  of  this  technique 
is  given  and  an  outline  given  of  the  design  methodology. 

1.  INTRODUCTION 

In  a  number  of  recent  papers  it  has  been  shown  that  H 
optimization  can  be  applied  to  the  problem  of  robust  stabiliza¬ 
tion  of  unstructured  uncertainty.  (See  Francis  [4],  and  the  ref¬ 
erences  therein.)  This  paper  describes  a  design  technique  based 
on  the  robust  stabilization  of  a  particular  representation  of  un¬ 
structured  uncertainty,  namely  that  of  normalised  ‘stable  fac¬ 
tor’  or  ‘coprime  factor’  perturbations.  (See  Vidyasagar  [15],  for 
a  suitable  introduction  to  coprime  factors,  and  [7]  for  a  com¬ 
plete  solution  to  the  normalized  coprime  factor  robust  stabiliza¬ 
tion  problem.)  The  design  technique,  which  also  allows  perfor¬ 
mance  objectives  to  be  incorporated,  has  two  stages:  (1)  A  loop 
•shaping  approach  ‘shapes’  the  nominal  plant  singular  values  to 
give  desired  open-loop  properties,  and  (2)  The  normalized  co¬ 
prime  factor  robust  stabilization  technique  mentioned  above  is 
used  to  stabilize  the  shaped  plant. 

2.  THE  NORMALIZED  LEFT  COPRIME 
FACTOR  ROBUSTNESS  PROBLEM 

We  firstly  summarize  the  main  results  of  [6],  [7]  for  the  nor¬ 
malized  LCF  robust  stabilization  problem.  It  has  been  shown 
[14],  [16]  that  an  attractive  way  of  representing  unstructured 
uncertainty  in  a  plant  is  via  coprime  factor  perturbations.  That 
is,  if  the  nominal  plant  is 

G  =  M~lN  (2.1) 

then  a  perturbed  plant  is  written 


Figure  2.1  Coprime  Factor  Uncertainty 

where  K  is  chosen  over  all  controllers  which  stabilize  G.  The 
solution  for  the  largest  achievable  c  (=  tmat)  is  generally  iter¬ 
ative  [16],  but,  if  the  left  coprime  factorization  (LCF)  of  G  is 
normalized ,  meaning 

+  N(ju>)N(ju>y  =  I  for  all  w  (2.4) 

then  it  is  possible  to  show  ([6],  [7])  that  a  maximum  value  of  e 
can  be  obtained  by  a  non-iterative  method,  and  is  given  by 

em„  =  (l-||[M,fV]||JH)0'5  (2.5) 

where  ||*||H  denotes  the  Hankel  Norm,  and  em0x  is  called  the 
maximum  liability  margin. 

Remark  2.1  It  can  be  simply  shown  that  the  problem  in  (2.3) 
is  equivalent  to  the  ‘four  block’  problem: 

K  (/  -  GK)-'\G,r) 
l  J  oo 

and  hence  this  also  has  a  minimum  solution  given  by  (2.5). 

3.  LOOP  SHAPING  METHODS 


GA=(M  +  AM)-'(A'  +  Ay)  (2.2) 

where  M,  N  is  a  left  coprime  factorization  (L.C.F.)  of  G,  and 
A/v  are  stable,  unknown  transfer  functions  representing  the 
uncertainty  and  satisfying  ||[Am,  A/v]^  <  e  where  e  >  0  (See 
Fig  21)  A  design  objective  is  then  to  find  a  feedback  controller 
A' which  stabilizes  all  such  G&  for  a  given  e.  Following  [16], 
this  can  be  rewritten  in  the  framework  of  an  H «,  optimization 
problem:  Find  a  stabilizing  conti  oiler  K  such  that 


In  feedback  design  many  performance  and  robust  stability 
objectives  can  be  written  as  requirements  on  the  maximum  sin¬ 
gular  values  of  particular  closed-loop  transfer  functions.  The 
principal  idea  of  ‘Loop  Shaping’  is  that  the  magnitude  (or  max¬ 
imum  singular  values)  of  these  closed-loop  transfer  functions  can 
be  directly  determined  (over  appropriate  frequency  ranges)  by 
the  singular  values  of  the  corresponding  open-loop  transfer  func¬ 
tion.  (The  reader  is  referred  to  [3]  for  a  comprehensive  introduc¬ 
tion  to  loop  shaping  methods.) 


I[  (I-GK)-'M-' 
\[K(I  -GK)-'M-' 


For  example,  for  a  plant,  G,  and  controller,  K ,  if  o{GK )  » 
; 2.3)  1  (typically  at  low  frequency)  then 
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4.  THE  DESIGN  PROCEDURE 


Figure  3.1  Loop  Shaping  Specifications 

?((/  -  GK )-' )  <  1  /a(GK)  (3.1) 

cf((J  -GK)'lG)  Zo{G)!g_{GK)  (3.2) 

and  if  qfGK)  C  1  (typically  at  high  frequency)  then 

o(K{I  -  GK )-> )  <  a(K)  (3.3) 

o{GK{I  -  GK)~l)  £  v(GK).  (3.4) 


For  good  performance  we  require  <?((/  -  GK)-1)  and 
a((l  —  GK)~'G)  to  be  small  (particularily  at  low  fre¬ 
quency)  and  for  good  robust  stability  properties  we  require 
o{K(I  -  GK)'1 )  and  a(GK(l  -  GK)-1 )  to  be  small  (at  high 
frequency  in  particular).  A  typical  closed-loop  design  specifica¬ 
tion  can  therefore  be  illustrated  as  in  Figure  3.1,  and  the  desired 
closed-loop  behaviour  can  be  achieved  by  manipulation  of  the 
open  loop  gains,  c(GK),  qfGK). 

However,  this  open-loop  shaping  approach  is  complicated  by 
the  need  to  ensure  stability  of  the  resulting  closed-loop  system. 
This  requires  that  plant  phase  properties  also  be  considered,  and 
the  loop  shape  can  be  shown  to  be  limited  by  such  stability 
requirements.  This  is  examined  in  (lj  for  the  SISO  case,  and  in 
[3]  for  a  MIMO  extension.  Further,  these  requirements  are  even 
more  restrictive  if  the  nominal  plant  has  RHP  poles  or  zeros. 
(See  [12]  for  example.) 

A  loop-shaping  approach  that  is  somewhat  simpler  from  the 
designers  point  of  view  is  that  used  in  the  Loop  Transfer  Re¬ 
covery  (LTR)  method  in  LQG  design.  (See  [9]  and  [2].)  In 
this  method,  the  designer  specifies  a  desired  singular  value  loop- 
shape,  and  the  guaranteed  stability  properties  of  the  LQG  com¬ 
pensator  ensure  stability.  LTR  however,  cannot  systematically 
deal  with  plants  with  RHP  zeros  (see  [13]),  and  is  limited  in  that 
it  can  only  guarantee  performance  and  robust  stability  proper¬ 
ties  at  either  plant  input  or  plant  output. 

The  design  technique  that  is  proposed  in  this  paper  is  sim¬ 
ilar  in  philosophy  to  LTR:  The  designer  specifies  a  desired  loop 
shape  and  then  the  'shaped'  plant  is  further  compensated  by  a 
controller  using  the  normalized  LCF  robust  stabilization  method 
of  Section  2  to  ensure  closed-loop  stability. 


4.1  Outline  of  the  Design  Procedure 

We  will  now  formally  state  the  design  procedure  that  was 
proposed  in  Section  1.  The  objective  of  this  approach  is  to  in¬ 
corporate  the  simple  performance/robustness  trade-off  obtained 
in  loop  shaping,  with  the  guaranteed  stability  properties  of  H ^ 
design  methods. 

‘The  Loop  Shaping  Design  Procedure  (LSDP)’ 

(1)  Loop  Shaping  -  using  a  precompensator ,  W1;  and/or  a  post- 

compensator,  W 2,  the  singular  values  of  the  nominal  plant 
are  ' shaped ’  to  give  a  desired  open-loop  shape.  The  nominal 
plant,  G,  and  ‘shaping  functions '  ,  Wj  are  combined  to 

form  the  ‘shaped  plant’,  Gs  where  Gs  =  WjGW, .  (See 
Figure  f.la.)  We  assume  that  YV j  and  Wj  are  such  that 
G s  contains  no  hidden  modes. 

(2)  Robust  Stabilization  -  a  feedback  controller,  KXl  which  ro¬ 
bustly  stabilizes  the  normalized  left  coprime  factorization  of 
Gs,  with  stability  margin  e  (=  J~x),  is  synthesised,  using 
the  approach  outlined  in  Section  2.  (See  Figure  f.lb.J 

(SJ  The  final  feedback  controller,  K,  is  then  constructed  by  com¬ 
bining  the  Hoo  controller,  Kx,  with  the  shaping  junctions,  W1 
and  W2  such  that  K  =  K^W] .  (Sec  Figure  4-lc.) 

Note  that,  in  contrast  to  the  classical  loop  shaping  approach, 
the  loop  shaping  here  is  done  without  explicit  regard  for  the 
nominal  plant  phase  information.  That  is,  closed-loop  stability 
requirements  are  disregarded  at  this  stage.  Also,  the  robust  sta¬ 
bilization  is  done  without  frequency  weighting.  The  parameter  e 
can  be  seen  as  an  indicator  of  the  success  of  the  loop  shaping, 
where  we  note  by  (2.5),  e  <  1  always.  A  small  value  of  c  (c  <  1) 


(a) 


<b> 


Figure  4.1  The  Design  Procedure  (a)  Shaped  Piaul,  (b) 
Robust  Stabilization,  (c)  Final  Controller. 
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in  Stage  i2l  always  indicates  incompatibility  between  the  speci¬ 
fied  loop  shape,  the  nominal  plant  phase,  and  robust  closed-loop 
stability  Hence  the  loop  shaping  of  Stage  { 1 )  is  connected  to 
the  robust  stabilization  of  Stage  (2)  via  the  indicator  < 

A  typical  design  works  as  follows:  The  designer  inspects  the 
open-loop  singular  values  of  the  nominal  plant,  and  shapes  these 
by  pre  and  'or  post  compensation  until  nominal  performance 
land  possibly  robust  stability)  specifications  are  met.  (Recall 
that  the  open-loop  shape  is  related  to  closed-loop  objectives.)  A 
feedback  controller,  Abo,  with  associated  stability  margin  (for 
the  shaped  plant)  <  <  fmai,  is  then  synthesized.  If  e  is  small, 
then  the  desired  performance  is  incompatible  with  robust  stabil¬ 
ity  requirements,  tend  the  loop  shape  should  be  adjusted  accord^ 
mgly,  and  Kx  re-evaluated.  (In  Section  5,  the  explicit  depen¬ 
dence  of  performance  and  robust  stability  objectives  on  <  will  be 

shown.  This  procedure  is  both  simple  and  systematic,  and  only 
assumes  knowledge  of  elementary  loop  shaping  principles  on  the 
part  of  the  designer. 

4.2  Comments  on  the  Design  Procedure 

The  following  comments  refer  to  the  design  procedure  given 
:::  Section  4  1. 

1.  In  the  Loop  5, .aping  Design  Procedure  we  have  interpreted 
’he  measure  >  as  a  design  indicator  rather  than  as  a  specific 
■  •.ability  margin  for  the  notional  perturbations  on  the  normalized 
LCF  of  the  shaped  plant:  In  Section  5  we  will  show  that  any 
-'abuizing  controller  achieving  (  <  1  will  lead  to  deteriorate u  in 
•he  achieved  loop-shape  com;  "  ’d  with  the  specified  loop  shape 
at  low  or  high  frequencies,  or  will  imply  poor  robust  stability 
properties  in  the  -elected  cross-over  frequency  region. 

2.  The  procedure  is  applicable  to  stable  or  unstable,  minimum 
or  non  minimum  phase  plants,  provided  they  satisfy  the  mini¬ 
mum  requirement  for  any  design  •  that  is,  there  are  no  unstable 
hidden  modes.  In  particular,  if  a  plant  is  non- minimum  phase, 
the  standard  performance  restrictions  will  still  exist  in  the  de¬ 
sign  procedure:  It  is  well  known,  (see  [8j  or  [5]  for  example) 
that  the  presence  of  a  non  minimum  phase  (RHP)  zero  in  the 
nominal  plant  limits  the  achievable  bandwidth  and  restricts  low 
'requeney  behaviour  In  the  Loop  Shaping  Design  Procedure,  a 
•  ■op  -hape  which  is  incompatible  with  these  generic  strictions 
will  lead  to  a  'small'  t  value  in  Stage  (2).  To  illustrate  this, 
- '.-is;  ier  the  following  example: 

Let  a  nominal  plant  be  given  by  G  =  ) ,  and  we  select 

U'  -  La  constant,  as  the  shaping  function.  The  shaped  plant  is 
•hen  Co  -  IV G  -  *'  ~ .  By  ;8!  we  know  that  the  closed-loop 

bandwidth  has  a  practical  limit  of  1  r/s  because  of  the  RHP  zero 
at  ,«  =  i  in  G.  Noting  that  the  crossover  frequency  of  the  open- 
’oop  shape  reflects  the  closed-loop  bandwidth,  we  now  select 
bur  values  of  k.  the  largest  yielding  a  crossover  frequency  that 
:s  incompatible  with  the  bandwidth  limitations.  The  resulting 
t  values  after  Stage  (2)  of  the  Loop  Shaping  Design  Procedure 
ar* 

1=1  <=06814 

k  =  1  <  -  0.5038 

*  c  .  0  2177 

k  =-  20  ■  --  2  0084 

Th'-  al ue  of  t  associa’ed  with  k  -  20  .  'vrtr'-nieiy  small. 
indica*.:r,g  an  ‘incompatible  loop  shape  1  This  is  confirmed  by 


io1  10*  ue  10'  to1  ice 

trtq  freq 

Specified 

Achieved 


Figure  4.2  Specified  and  Achieved  Loop  Shapes-RHP 
Zeros  Example 

comparing  the  specified  and  achieved  loop  shapes  for  each  case 
in  Figure  4.2. 

In  the  next  section,  we  give  theoretical  justification  for  using 
the  loop  shaping  design  procedure 

5.  PROPERTIES  OF  THE  DESIGN  METHOD 

5.2  Controller  Magnitude  Bounds 

In  the  previous  section  we  specified  the  desired  loop  shape 
by  WjGWi  (Fig.  4.1a),  but  the  actual  loop  shape  achieved  is  in 
fact  given  by  W1  A’oo  WG  at  plant  input  (point  (i),  Fig  4.1c), 
and  G'V\h'x>W2  at  plant  output  (point  (ii),  Fig.  4.1c).  We 
will  now  show  that  the  degredatinn  in  the  loop  shape  caused  uy 
the  Hx  controller  Kx  is  limited  at  those  frequencies  where  the 
desired  loop  shape  is  sufficiently  large  or  sufficiently  small: 

5.1.1  Low  Frequency  Behaviour 

At  lowfiequency  (in  particulars  €  (0,w;))  the  deterioration 
in  loop  shape  at  plant  output  cam  be  obtained  by  comparing 
<t| GW,KX,W2)  with  riW'jGW,!.  Note  that: 

a(GK)  =  a  {GW,  K^W2)  y  u(  W2GW,  )o(A'00)/c(  W,  ) 

(5.1) 

Sirinlarily,  for  loop  shape  deterioration  at  plant  input  we  compare 
ct.(W,K'kWjG)  with  ^(W'jGWj)  and  we  have 

o(KG)  =  r'W,  Abo  W2G)  >  o{W%GWx)o(Kx)lciWi) 

(5.2) 

In  each  case,  o(  A'.*,)  is  required  to  obtain  a  hound  on  the  dete¬ 
rioration  in  the  loop  shape  at  iow  frequency.  ‘  ote  that  the  con 
dition  numbers  cl  W  !  and  c(  VVL )  arc  selected  by  the  designer, 
and  are  commonly  ol  order  one 

The  following  result  shows  that  nj  Kx  )  is  explicitly  bounded 
by  functions  of  >  and  rrtG’s),  the  minimum  singular  value  of  the 
shaped  plant.  We  a-ssume  here  that  the  shaped  plant  Gs  has 
an  equal  number  of  input--  and  outputs  For  convenience,  the 
parameter,  -  defined  t  C  '  1  will  be  u-.-d  in  the  following 
analysis 
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Theorem  5.1  Any  controller ,  Ka c.  satisfying 

jS  [  \  -GsK^y'M-'  <7 

where  i,V#l  A/,)  is  a  normalized  LCF  of  Gs,  nlso  satisfies 
aV'sO-))  (73-l) 


\Hl  -  1(1  +q2{Gs(joJ)))  r  7  1^[g7u^) 


-6.3) 


/or  d/.1  _■  juc/i  f/iaf 

0,GsO~-))  >  \A2  -  1. 


/or  a/i  Jtic/i  i/ict 

o(Gs{jui))  < 


n/t2”-  1 


Again,  if  w«-  consider  those  frequencies  at  which  ct(C»s) 

-  L  =  th<*n  we  have  the  following  result: 

vAJ-> 

Corollary  5.4  Following  the  notation  of  Theorem  5.3,  if 
a(G s(ju))  ^  f  then 


o(K^(j u/))  <  '/* 


(Lengthy  proofs  axe  omitted  in  this  paper,  and  the  reader 
is  referred  to  [1 1  j  for  hill  details.) 

The  main  implication  of  Theorem  5.1  is  that  the  bound  on 
'j:  h\c  j  depends  only  on  the  selected  loop  shape,  and  the  sta¬ 
bility  margin.  The  value  of  7  directly  determines  the  frequency 
range  over  which  this  result  is  valid  -  a  small  7  is  desirable,  as 
we  would  expect.  Further,  if  we  consider  those  frequencies  where 
'7i  Gs )  \/72  —  1,  we  have  the  following  asymptotic  result: 

Corollary  5.2  Following  the  notation  of  Theorem  5.1,  if 
l V  G s\J-^ ) )  VV  1»  then 

oSK^jui))  >  -4==. 

VT2  -  1 

Proof:  This  follows  immediately  from  equation  (5.3),  by  noting 
that  Z\Gj{j~-))  -  1  implies  that  S.2{Gs(ju))  3>  1,  and 

also  that  zjl  G ))t~r‘  -  1)1/2  >  7 2<l(Gs(;w)).  ■ 


5.1.2  High  Frequency  Behaviour 

At  high  frequency  (in  particulars  t  (w„,oo))  the  deterio¬ 
ration  in  plant  output  loop  shape  can  be  obtained  by  comparing 
aiGWx  K vjW?)  with  o(WjGWi).  Note  that,  analogously  to 
(5  i)  and  (5.2)  we  have 

oyGK)  =-  o'GW.K^W,)  <  a(W2GWiy<f(K^)c(W2).  (5.4) 

Similariiy.  the  correspond!,  g  deterioration  in  plant  input 
loop  shape  is  obtained  by  comparing  a(  Wj  W2  G)  with 
r'  VV; Cl IV;  j  where 


'  KG,  ,  a,  VV,  K.,W,G)  <  <f(  )a(K^)c(  W, ).  (5.5) 


Hence  in  each  case,  of  A'.*,)  is  required  to  obtain  a  bound  on 
the  deterioration  in  the  loop  shape  at  high  frequency.  In  an 
analogous  manner  to  Theorem  5.1,  we  now  show  that  cffATo)  is 
explicitly  hounded  by  functions  of  7,  and  o(Gs ),  the  maximum 
singular  value  of  the  shaped  plant. 

Theorem  5.3  Any  controller,  K^,,  satisfying 


r  A'-x 
I  / 


1  ■'  / 
J 


GsK^r'M; 


■30 


£  7 


where  (AT.  M.,  a  normalized  LCF  of  Gs,  also  satisfies 


Proof:  This  follows  immediately  from  (5  6)  by  rioting  that 

alGs)  <  -v"1; implies  cr'(Gs )  <  1,  and  -.;r|G's)  -K  17*  - 
Va’-i 

I)1'*.  ■ 


Remark  5.5  The  approximate  results  in  C.toII.uks  5  2  .aid 
5.4  show  that  at  frequencies  where  q_[Gs)  S>  1  or  <jtGsl  «u  !. 
the  deterioration  in  the  loop  shape  due  to  A  is  hounded  by 
a  function  of  7  (or  equivalently  e)  only.  Noting  that  7  'small' 
(implying  (  'large')  indicates  a  minimal  deterioration  in  the  loop 
shape,  it  is  confirmed  that  e  indicates  the  compatibility  between 
the  specified  loop  shape  and  closed-loop  stability  requirements 

Remark  5.6  Note  that  in  Theorem  5.3  it  is  not  necessary  to 
assume  that  Gs  is  a  square  transfer  function  matrix  Hence, 
high  frequency  loop  shape  properties  ran  be  guaranteed  for  a 
plant  of  any  dimension. 

We  have  shown  in  this  section  that  the  values  of  7  (alter 
natively  ()  achieved  in  the  Loop  Shaping  Design  Procedure  will 
directly  affect  the  singular  values  of  the  //„  controller  AT*,  Al¬ 
though  the  precise  relationship  between  7  and  the  shaping  func¬ 
tions  4V [  and  VV2  is  not  known,  we  have  shown  that  a  large  7 
value  indicates  incompatibility  between  the  selected  loop  shape 
and  the  closed- loop  stability  reqiurements.  Such  an  incompat 
ibility  can  lead  to  extensive  deterioration  in  the  specified  loop 
shape.  Under  such  circumstances,  the  designer  would  be  re 
quired  to  ‘relax’  the  original  loop  shape  specification,  unci,  a 
more  compatible  loop  shape  in  achieved 

5.2  Using  7/f  as  a  Design  Indicator 

In  Section  5.1  we  evaluated  bounds  on  the  ma.-imiirii  loop 
shape  deterioration  at  low  and  high  loop  gains,  for  a  given  -j  — 
e~'),  and  a  specified  loop  shape.  In  this  section  we  lake  a  ui 
trary  approach:  Given  any  stabilizing  controller.  AT  .  there  exists 
a  frequency  such  that 


(/  -  GsK.r 


0- o' 


and  we  evaluate  the  minimum  deterioration  in  loop  shape  that 
can  be  expected,  if  any,  for  this  7,  and  examine  the  effect  o i  7 

on  robust  stability  at  this  frequency.  VVe  sumii.ari/e  lie  .<  j.i,  ,L 
in  the  following  theorem 

Theorem  5.7  Let  ',-uj  be  the  optimal 

for  the  normalized  LCF  robust  stabilization  pt.ihit  si.  r. 

(2.5)  Then,  for  any  stabilizing  tontrollrr  NT.  t -■  •  .is. 
frequency,  ^z0,  surh  that 


G ,,  A0 ! 


Af. 


n 


and  hence 


(\)  There  ezists  a  perturbed  system 


G±  =  Wt-l(M,  +  +AN)Wr‘ 


with 

!![A,v.AM]IL  (5.9) 

which  destabilizes  the  system.  Further , 
i  n)  if  >  r:i,n  -  t/ien 


V>m,n  -  U1  +  g‘(GsQw,,)))  +  rlmg  Gs(jUo)) 
(7^in-l)22(Gs(j^))-l 


ana 

fliij  >/ j'Gsljw',))  <  (7m,„  -  1)1/2,  then 


7m, n  ~  1  ~  <j2(Gs\]~o  l) _ 

1(1  +  a2(Gs{jui0)))  +  ?^,„aiGs(juj0)) 


Remark  5.8  Theorem  5.7  shows  that  •>/«  is  a  design  indicator 
for  the  entire  frequency  range:  Firstly,  if  the  coprime  factor  ro¬ 
bust  stability  is  poor  (that  is  t  <  1  or  7  >  1)  at  a  frequency 
where  N,  and  M,  are  of  comparable  size,  then  (•'  shows  that 
only  a  smali  relative  perturbation  on  either  TV,  or  iVf,  is  per¬ 
mitted.  However,  if  M,  ~S>  TV,,  then  the  notional  destabilizing 
perturbation  could  correspond  to  a  very  large  relative  pertur¬ 
bation  in  .V, ,  implying  that  robust  stability  properties  may  in 
fact  be  acceptable  despite  7  >  1.  However,  in  this  case,  we 
have  ?(G,-j  <  1  and  part  (iii)  of  Theorem  5.7  applies,  show¬ 
ing  that  the  loop  gain  is  necessarily  and  undesirably  increased. 
Conversely,  if  M,  <  TV,,  then  the  notional  destabilizing  pertur¬ 
bation  could  correspond  to  a  very  large  relative  perturbation  in 
M,.  again  implying  that  robust  stability  properties  may  in  fact 
be  acceptable  despite  7  >  1.  However,  in  this  case,  we  have 
niGs  i  3>  1  and  part  (ii)  of  Theorem  5.7  applies,  showing  that 
the  loop  gain  is  necessarily  and  undesirably  decreased. 

In  the  next  section  we  examine  the  closed-loop  behaviour 
achieved  using  the  Loop  Shaping  Design  Procedure,  and  confirm 
that  the  loop  shaping  approach  gives  guarantees  on  particular 
cioced-loop  objectives. 


5.3  Behaviour  of  Standard  Closed-Loop  Objectives 

In  Section  3  we  stated  that  a  feature  of  the  classical  loop 
shaoing  design  approach  is  that  it  is  possible,  by  open-loop  sin¬ 
gular  value  shaping,  to  ensure  that  a  number  of  standard  closed- 
loop  design  objectives  are  'well  behaved’.  The  following  result 
demonstrates  that  this  is  also  the  case  for  the  Loop  Shaping 
Design  Procedure  outlined  in  Section  4.1: 

Theorem  5.9  Lei  G  be  the  nominal  plant  and  let  K  = 
VV)  K~c  \V2  be  the  associated  controller  obtained  from  the  Loop 
Shaping  Design  Procedure  of  Section  4.1.  Then 


and 


0 (K  I-  GK)~l)<  -,a<M,WWt)a(W2) 
o\(I  -  GK)~X )  <  7<f(  M,  )c(  W2 ) 
a~  K  1  -  GK  r'G)  <  7<f( N, )c(  W, ) 


ffh  I  -  GK)~ 


G)  < 


7  er{  IV, )  _ 
at  Wt  W2 ) 


(5.10) 

(5-il) 

(5.12) 

(5.13) 


a((I-  KG)-')  <  l+io(N,)c(Wi)  (5.14) 

c(G(I-KG)~lK)<  1  +  7  <j(M,)c(W2)  (5.15) 


where  (IV,, IW,)  is  a  normalized  LCF  of  Gs  =  W'jGW'i,  and 
c(»)  denotes  the  condition  number. 

Remark  5.10  Note  that  because  (IV,,  M,)  is  a  normalized  LCF 
of  Gs  we  have  that  o( IV, )  <  1  and  a(M,)  <  1  for  all  frequen¬ 
cies.  Noting  also  that  7  has  a  finite  value  ( typically  1  <  7  <  5  in 
practice),  and  that  the  shaping  functions  are  selected  by  the 
designer,  then  it  can  be  seen  that,  by  (5.10)  -  (5.15),  all  of 
the  closed-loop  objectives  are  guaranteed  to  have  bounded  mag¬ 
nitude.  We  say  that  in  this  case  the  objectives  are  well  be¬ 
haved.  For  the  simple  case  of  =  J  and  W2  =  L  it  can 
be  readily  seen  that  in  (5.10)  -  (5.13),  IV,  (=  TV)  and  M, 
(=  M )  provide  frequency  weighting.  The  frequency  shapes 
are  ‘natural’  for  many  problems:  for  example  at  frequencies 
where  o(M,)  is  small  (indicating  a  nominal  plant  pole  near  the 
imaginary  axis),  then  a((I  -  GK)~' )  and  a(K(I  -GK)~l) 
are  small.  So  the  behaviour  of  the  standard  closed-loop  transfer 
functions  (I  -  GK)~' ,  K(I  -  GK)~l ,  K(I  -  GK)~lG,  and 
(I-GK)-'G  is  compatible  w-ith  sensible  closed- loop  design. 

5.4  Bounds  on  the  Normalized  Coprime  Factors 


We  now  state  a  technical  result  which  demonstrates  that 
c(N,)  and  a(M,)  are  related  to  the  nominal  plant  G  and  the 
shaping  functions  W\  and  \V2 . 


Lemma  5.11  Let  the  shaped  plant,  Gs  = 
normalized  LCF  given  by  Then 

■■  W2GxV j,  have  a 

o(N3)  = 

/  aHW2  gw,)  ' 

iG 

)  (5.16) 

^l  +  ^WjGWi), 

and 

V(M.)  =  1 

f  1 

1/2 

)  •  (5.17) 

U  +q}(W2GW,)J 

Remark  5.12  The  value  of  this  Lemma  is  that  the  bounds 
on  the  closed-loop  objectives  in  Theorem  5.9  can  be  rewritten 
in  terms  of  7,  G,  Wj,  and  Wj  only.  It  can  now  be  clearly 
seen  bow  the  loop  shaping  influences  the  closed-loop  properties. 
Noting  that  if  <7(VV2GWj)  >  1  then  a(M,)  nr  l/e(W,2GW'1), 
<t(IV,)  ~  1  and  that  if  a (WoGW^)  <  1  then  <f(jVf, )  ~  1 
a(N,)  nr  ct( VV^GWi ),  then  the  bounds  in  Theorem  5.9  can  be 
re-evaluated  to  show  the  effects  of  W\  and  W2  on  closed-loop 
behaviour  in  frequency  regions  of  high  and  low  loop  gain. 


6.  SUMMARY  OF  RESULTS 

In  this  paper  we  have  incorporated  the  normalized  LCF 
robust  stabilization  problem  into  a  loop  shaping  based  design 
technique.  This  enables  both  performance  and  robust  stability 
objectives  to  be  traded-off,  and  preserves  the  exact  solution  as¬ 
sociated  with  this  particular  Hoc  pioblem.  The  design  method  is 
straightforward  and  systematic,  incorporating  only  the  simpler 
aspects  of  ‘classical’  loop  shaping.  The  normalized  LCF  robust 
stabilization  problem  has  been  shown  to  be  particularily  well 
suited  for  this  approach  because  the  Hoc  controller  synthesized 
causes  only  a  limited  deterioration  of  the  specified  loop  shape 
Coprime  factor  model  reduction  techniques  ran  be  simply  incor¬ 
porated  into  this  framework  As  a  final  comment,  it  should  be 
noted  that  the  Loop  Shaping  Design  Procedure  is  restricted  to  a 


78 


particular  set  of  Hx  objectives.  However,  it  can  be  shown  that 
the  design  technique  yields  simple  and  effective  controllers.  (See 

[10]  and  [11].) 
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Abstract 

The  problem  of  maintaining  the  Hoc  norm  of  a  standard  closed  loop  below  a 
prespeafied  tolerance  level  whilst  maximizing  an  entropy  integral  (at  a  point 
in  the  right  half  plane)  is  posed  and  solved  by  way  of  the  equivalent  error 
system  distance  problem.  All  error  systems  with  infinity  norm  below  the 
tolerance  level  are  parametrized  as  the  linear  fractional  map  of  an  all  pass 
matrix  and  an  arbitrary  stable  contraction.  We  derive  the  maximum  entropy 
choice  for  the  contraction  and  a  value  for  the  maximum  entropy.  For  the 
maximum  entropy  at  infinity,  it  is  proved  that  the  arbitrary  contraction  must 
be  set  to  zero  (the  ‘central’  solution);  the  maximum  entropy  in  this  case  is 
an  explicit  formula  in  terms  of  the  realization  of  the  error  system.  Some 
motivational  remarks  are  made  and  links  between  entropy,  H2  norms  and  H3 
optimal  control  are  giv-n. 


Notation 

'  ^  The  open  right  half  plane. 

7Z  (Prefix)  Real -rational- 

KHj  Hardy  space  of  real-rational  transfer  function 

matrices  analytic  and  square-integrable  on 
vertical  lines  in  the  right-half  plane. 

TJHoo  Hardy  space  of  real-rational  transfer  function 

matrices  analytic  and  bounded  in 
the  right  half  plane. 

G*  :=  G(-$)T,  the  parahermitian  conjugate. 

<7-,(G)  :=  \1/J{G’G),  the  Xth  singular  value  of  G. 

iiG'::,  =  !rr rrace(G'0'w)G0'u-))<iv]  ,  H,-nor m. 

!!G||„  .=  sup„ff„„(G(;w)),  H<«, -norm. 

X  =  Ric  (  £.c  -BAB.  )  '"is  the  stabilising  solution  of  the  Riccati  equation 
XA  *  A’X  -  XBB’X  +  C'C  =  0 


1  We 


l  )h3 


MS) 


Denotes  maximum  entropy 
Denotes  Ho,  optimal. 

Denotes  H,  optimal. 

:=  D  +  C(jI  -  A)_,B,  state-space  realization. 


A  square  transfer  function  matrix  G(s)  is  said  to  be  all  pass  if 


Figure  1.  The  standard  configuration. 

The  closed  loop  transfer  function  matrix  G  from  w  to  z  is 
given  by  the  linear  fractional  map 

G  =  ^(P,K) 

:=P„  +P,aK(I  — P„K)--P„  (1) 

of  the  appropriately  partitioned  standard  plant  P  and  the  con¬ 
troller  K.  The  optimal  H„  control  problem  is  to  find  a  stabi¬ 
lizing  controller  Kop,  which  minimizes  the  infinity-norm  of  this 
transfer  function  i.e.,  K„p,  satisfies 

«n/{||^(P,K)|L  :  K  stabilizes  P}  :=  7,„ 

=  |Pr(P,K.„)||00.  (2) 

Motivated  by  the  belief  that  optimal  H„,  control  is  not  always 
appropriate  we  consider  here  the  suboptimal  problem  obtained 
by  relaxing  the  infimum  in  (2)  to  the  bound 

ll^(P,K)|L<7  (3) 


G*(ju)G(;w)  =  I,  Vw. 

The  Laplace  transform  variable  s  will  be  suppressed  to  keep  typography  sim¬ 
ple.  Dependence  of  transfer  functions  on  •  is  understood. 

I.  Introduction 

This  paper  is  concerned  with  suboptimal  control  with 
a  maximum  entropy  criterion.  The  theory  of  optimal  H*,  con¬ 
trol  has  received  much  attention  over  recent  years;  for  full  de¬ 
tails  the  interested  reader  is  referred  to  Francis  [6]  and  the  refer¬ 
ences  therein.  Figure  1  illustrates  the  usual  configuration,  where 
the  ‘standard  plant'  P,  consisting  of  the  actual  plant,  suitable 
weighting  functions  and  interconnections,  maps  exogenous  in¬ 
puts  w  and  control  inputs  u  to  controlled  outputs  z  and  mea¬ 
sured  outputs  y.  As  is  usual  in  H.  control  problems,  we  assume 
'-hat 


P  = 


Sit 
S3  I 


( 


Pl. 

P>. 


m,  >  Pi ,  Pi  >  m, . 


where  7  >  7,,,.  In  general,  there  is  a  class  of  controllers  which 
satisfy  this  bound;  such  nonuniqueness  is  dealt  with  in  this  paper 
by  specifying  that  the  entropy  of  the  closed  loop  transfer  function 
must  be  maximized. 


The  entropy  is  defined  as  follows.  For  any  transfer  function 
matrix  G  which  satisfies  ||G|j„  <  7  the  entropy  of  G  at  a  point 
J0  6  C+  is  defined  by 


I(G;7;Jo) 


r  r  ln)det(l  —  7_1  G"  (jw)G(jui))| 

J-ao 


Re  s„ 
l«o  -  M 


duj . 


w 

This  definition  is  equivalent  to  that  of  Arov  and  Krein  [  1 ) ,  [2] 
except  for  an  extra  Re  s0  term  to  ensure  non  zero  entropy  when 
s„  — >  00.  It  is  easily  seen  that  the  entropy  is  well-defined  (since 
IIGHc  <  7  implies  0  <  I  —  7" 1 G*  (jui)G{jai)  <  I  )  and  non 
positive;  that  I(G;  7;  s„ )  =  0  if  and  only  if  G  =  0;  and  that  the 
entropy  of  G  is  invariant  under  unitary  scaling  of  G 


RPMOco  1.7/00  eonr.oocot  -  -  .ryoc  icer 


HO 


Our  maximum  entropy  control  problem  is  then: 
hind,  "nt  nf  all  rontmllrrs  K  which  stabilize  I*  and  satisfy 

1 1 -HP  K )  ||  ,  <  7  (5) 

the  K  which  maximizes  the  entropy  I(G;  7,'  30)  of  the  closed  loop 
transfer  function  matrix  G  =  .F(P,  K)  at  a  point  J,  6  C,  . 

Maximum  entropy  has  been  studied  in  a  wide  variety  of 
contexts;  its  applications  to  extension  problems  (Arov  and  Krein 
1],  [2j)  and  to  contractive  interpolants  (Dym  and  Gohberg  [4j, 
[5])  are  pertinent  here  -  we  use  an  adaptation  of  the  method 
of  [2].  The  use  of  maximum  entropy  in  Hoc  control  has  been 
considered  by  Limebeer  and  Hung  [9]  for  the  ‘one-block  case' 
where  bo'h  P,>  and  P,,  are  square.  No  v_ch  assumptions  are 
made  here,  the  analysis  is  valid  for  ‘one-,’  ‘two-,’  and  ‘four- 
block’  problems. 

As  with  optimal  H«,  problems  (see  [6])  we  approach  the 
problem  by  reducing  our  original  problem  to  a  ‘distance  prob¬ 
lem'  .  To  do  this,  we  use  the  parametrization  of  all  stabilizing 
controllers  of  Youla  et  at.  [I2j  to  reduce  (5)  to  the  equivalent 
model-matching  problem  of  finding  Q  €  such  that 

||T,  +  T,QTj  Ha,  <  7  (6) 

and  then  exploit  the  unitary  invariance  of  the  H^,  norm  to  re¬ 
duce  (6)  to  the  distance  problem 

S(£:  ^kh-  <7) 


where  R  =  (  j*'1 1  ?'1  J  |  is  antistable  and  is  known  in  terms  of 

\  J  3  / 

the  standard  plant  P. 


If  we  define  the  error  system  E  by 


p  /  ^n  fin  \ 
_  \ Rsi  R,3  +  Q ) 


we  know  that  I(G;7l«o)  =  I(E;  -y;  )  because  entropy  is  unitar  - 

ily  invariant,  therefore  the  closed  loop  transfer  function  G  and 
the  error  system  E  have  the  same  entropy,  allowing  us  to  solve 
our  original  closed  loop  problem  (5)  by  solving  the  following  error 
system  distance  problem: 

Find  the  error  system  E  as  defined  m  equation  (8)  which  satisfies 
||E|U  S  7  ond  maximizes  I(E;  7;s0 ). 

The  arrangement  of  the  paper  is  as  follows.  In  the  next 
section  we  briefly  motivate  control  and  maximum  entropy. 
Section  III  contains  the  main  results.  By  parametrizing  all  solu¬ 
tions  of  the  error  system  distance  problem  we  are  able  to  derive 
the  unique  maximum  entropy  solution  together  with  a  value  for 
■  t.-  ntropj  This  is  firstly  done  for  tin  general  case  of  entropy  at 
any  s„  6  C,  and  then  for  the  important  special  case  of  s„  — >  00. 
This  latter  case  yields  particularly  explicit  and  appealing  results. 
Proofs  of  certain  lemmas  have  been  relegated  to  the  appendix 
whenever  the  inclusion  of  the  proof  in  the  main  body  of  the  text 
would  be  int.  J3ive. 

II.  Motivation 

Here  we  briefly  state  some  relevant  background  details.  Re¬ 
call  that  we  want  our  controller  to  stabilize  P  and  keep  the 
infinity-norm  of  the  closed  loop  F(P,K)  below  a  level  7  (where 
7  >  7«pi  )•  Such  control  problems  lead  to  a  class  of  possible  con¬ 
trollers  and  so  there  is  scope  for  another  criterion.  Our  approach 
is  to  use  the  controller  which  maximizes  the  closed  loop  entropy. 
The  closed  loop  entropy  (4)  is  a  useful  measure  of  how  close 
G  =  TiP.K)  is  to  the  worst  case  of  G*(juj)G (jui)  =  7’ I  Vu;, 
where  the  bound  is  achieved  at  every  frequency,  for  in  this  case 
the  entropy  is  -  oc 


As  G  moves  away  from  this  worst  case  towards  G  -  0 
the  entropy  becomes  finite  and  decreases  in  magnitude  until  at 
G  -  0  the  entropy  equals  zero  Maximizing  the  entropy  at 
any  particular  7  is  an  effective  way  of  driving  G  away  from  the 
G"  [ju)G(juz)  =  71 1,  Vui,  case  towards  the  more  desirable  G  =  0 
case.  Of  course,  our  standard  plant  must  be  stabilized  as  well. 
Furthermore,  if  we  rewrite  the  entropy  as 

I(G;7;So)  =  i  r  X^'nll-7-V,1(G(;W))|[r^-|>  du, 
J-  00  ,  IK  -  j-j  I J 

(9) 

it  is  clear  that  ail  the  singular  values  <r,{G)  of  G  are  included, 
unlike  the  infinity  norm  which  depends  only  on  the  largest  sin¬ 
gular  value. 

The  term  [(Res0)/|s0  —  yui|)’  in  the  entropy  integral  is  a 
frequency  weighting  with  a  shape  dependent  on  the  position  of 
the  point  s0  in  the  right  half  plane.  In  order  to  obtain  real- 
rational  controllers  s0  should  be  a  real  numbe1-;  allowing  .s. 
makes  the  frequency  weighting  equal  to  unity  for  all  frequencies, 
a  notable  special  case  we  will  return  to  later. 

An  interesting  link  with  the  H3  norm  is  provided  by  the 
following  lemma. 

Lemma  2.1.  For  any  G  which  satisfies  ]| G j| „  <  7  we  have 
(,)  {-I(G;7,s0)}‘/j  >  ||G(fie  s0)/(s,  +  s)||, 

(it)  { — 1(  G ;  7;  oo)}’/J  >  ||G||j  if  G  is  strictly  proper 
and  equality  m  both  cases  is  achieved  when  7  ->  oc 
Proof:  Appendix.  AAA 

Part  (i)  of  this  lemma  shows  us  that  the  square  root  of  the 
magnitude  of  the  entropy  at  s0  is  an  upper  bound  on  a  frequency 
weighted  H,  norm,  whilst  part  (ii)  is  a  similar  result  for  the  en 
tropy  at  infinity  and  the  usual  (unweighted)  H,  norm,  if  it  exists 
In  both  (i)  and  (ii),  relaxing  the  H„  norm  constraint  entirely 
by  allowing  7  — *  00  gives  equality.  In  other  words,  the  frequency 
weighted  (respectively  unweighted)  H,  norm  minimization  prob¬ 
lem  is  exactly  the  maximum  entropy  at  s„  (resp.  at  do)  prob¬ 
lem,  with  7  — >  co.  Thus  the  imposition  of  the  maximum  entropy 
constraint  allows  us  to  use  7  to  trade  off  between  H„  optimal 
(7  — *  7,,i)  and  H,  optimal  (7  -v  00)  solutions. 

III.  Derivation  of  the  maximum 
entropy  solution 

In  this  section  we  solve  the  following  maximum  entropy  dis¬ 
tance  problem,  as  posed  earlier. 

Let  s0  e  C,  and  ~ 


R  _  Fl-">ll  |  Ru] 

Pit  \  1^11  / 


be  given,  where 


R’  £  ,  m,  >  pi ,  p,  > 


Define  the  error  system  E  by 


and  let 


’•(&:  r,R”q) 


7.,,  =  mf  { 1 1 E j ; ^  :  Q  t  7CH„ 


Then  for  7  >  7.,, ,  find  Q  6  7JH„  such  that  IJE!  , 
entropy  I(E;7;  sQ  )  15  maximized 
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We  have  seen  that  this  problem  is  equivalent  to  finding  a 
stabilizing  controller  K  which  keeps  ||^(P,  K)||.  <  7  and  max¬ 
imizes  the  closed  loop  entropy. 

Solution  proceeds  by  firstly  parametrizing  all  possible  E 
which  satisfy  the  bound  HE!!,*,  <  7.  Ball  and  Cohen  (3]  pro¬ 
vide  such  a  parametrization  in  terms  of  a  linear  fractional  map 
of  a  J-um tary  matrix  and  an  arbitrary,  stable  contraction  $ 
(that  is,  *  £  KH*  and  ||*||.  <  1),  but  it  is  more  convenient 
to  use  the  parametrization  of  Glover  and  Doyle  [8]  in  terms  of 
the  linear  fractional  map  of  an  all  pass  matrix  and  an  arbitrary, 
stable  contraction  <$.  By  adapting  the  method  of  [2]  we  are  able 
to  derive  the  unique  choice  of  *  which  maximizes  the  entropy 
and  a  value  for  the  maximum  entropy,  for  both  the  general  case 
of  s0  6  C+  and  when  j0  — *  00. 

III.l  The  general  case 

Here  we  solve  the  maximum  entropy  distance  problem  for 
arbitrary  s0  £  C+  and  E  proper,  but  not  necessarily  strictly 
proper.  The  class  of  error  systems  E  over  which  the  entropy 
must  be  maximized  is  parametrized  in  the  following  lemma. 
Lemma  3.1.  [8]  All  solutions 


E  = 


"M  -P2  ^ 

n-ntl  f  R-11 

">it  l  ^>i  +  Q 


with 

R‘  ■  Q  e  ftHoo  ,  m,  >  p3 ,  p,  >  m, , 

to  the  distance  problem  ||E||.  <  7  where  7  >  7.,,,  air  given  by: 


E  =  7^(R..  +  q..,n 


(13) 


Lemma  3.2.  Suppose  J  =  (j’j  j”)  is  all  pass,  *  £  7?H.  , 
||\&!U  <  1  and  det(l  -  J33*)  is  a  unit  in  72H.  .  Then 

I(7^(J,*);7;s0)  =  7T(*;l;So)-i  73I(J. . ;  1;  J.) 

-7J(*e  s0)/n|de!(I  —  (*. )*(«•  ))|.  (16) 


Proof:  Throughout  this  proof  take  s  =  jcj.  The  assumption 
that  det{ I  -  J334>)  is  a  unit  in  V.li„  together  with  the  all  pass 
nature  of  J  ensures  that  (Redheffer  [10]) 

(7^(J,*)]'[7^(J,*)]  <7*1.  (17) 

so  the  entropy  1(7^(J,  *);  7;  s0 )  is  well  defined.  Using 

J’(ju)JU^)  =  I  (18) 

in  block-partitioned  form,  it  is  straightforward  to  show  that 

I-7-’l7JF(J,*)]-  h^(J,*)] 

=  J3\  (I -J„ *]-•{!-***]  (I-J,,*)-1  j„.(19) 

From  this,  and  the  fact  that  for  any  square  real-rational  transfer 
function  matrix  G 

ln\det(G‘  G)|  =  ln|det(G)|  +  /n|de<(G‘  )|  =  2  /n|det(G)|,  (20) 

we  obtain 

/n|det(I  — 7-’[7^(J,*)n7^(J,*)l)l 

=  /n|det(I-  *-*)|  +  /n|det(j;iJ31)|  -2  fn|det(I  -  J„«)l- 

(21) 

Substituting  this  into  the  definition  of  I{7^(J,  *);7;  s0)  and  us¬ 
ing  the  (1,1)  block  of  (18)  to  write  J  j ,  J, ,  =  I  -  J*t  J, , ,  Vs  =  ju, 
it  follows  that 


where 

*  :  = 

Also, 

R-aa  -1"  Q„ 


"1  -Pa  2 

I(7^(  j ,  * ) ;  7;  ^0 )  =  7J  I(  *;  1 ;  )  +  7s  I(  j. . ;  1 ;  s„ ) 

(  0  *) 

,  *€71H„,  ||*|1.  <  1. 

/  fn|def(I  -  Jj3(jw)*(jw)| 

’  Res ,  l3 

.(•So  -  jut  1. 

_  nl  (  [R.a  +  Q.» 

)i  1  [R-a •  +  Qaa  jl  3 

) 

(14) 

"'ll  \  [R»«  +  Q.a 

|a  1  [R-ao  +  Qaa]aa 

) 

m,  ~p? 

Pa 

pi  -"*? 

”*3 

P\  -  *» al 

V'R,, 

7“  1  Ri  a 

R.. 

0  N 

—  m  1 1 

7” '  Rji 

7“‘  (Raa  +  Qaa) 

Roj  +  Qjj 

Q,. 

"»1  -Pll 

R„ 

R^a  +  Qaa 

Rjj  +  Qj> 

Q,. 

Pi ! 

V  0 

Q*  a 

Qn 

Q ../ 

du. 

(22) 

By  assumption,  det(I  —  J33*)  is  a  unit  in  72.H.  ,  which  per¬ 
mits  the  use  of  Poisson’s  Integral  Theorem  (Rudin  [11],  p343)  to 
evaluate  the  integral  in  (22),  giving 

I(7^(J,  *);  7;  a. )  =  7’  I(*;  1;  *. )  +  7*  I(  J. , ;  1;  *. ) 

-  7 ’(Re  S0  )ln\det(l  -  J33(s0  )*(s0  ))|  (23) 


as  claimed. 


AAA 


(15) 


Further,  ,Q,,  €  R.H.  ,  R,,  +  Q,.  is  ail  pass,  Q,»(oo)  =  0 
and  HQ,.!!.  <  1. 

State  space  realizations  of  R,y  and  Q,y  are  available  in  [8j, 
in  terms  of  the  realization  of  {*"  S'15 )  and  the  solutions  to  two 

'Rj|  Raa  t 

algebraic  Riccati  eq^Rions.  In  this  section  we  will  not  need  these 
realizations. 

The  next  lemma  relates  the  entropy  of  the  linear  fractional 
map  of  an  all  pass  matrix  J  and  an  arbitrary  stable  contraction 
*  to  the  entropy  of  ¥  itself. 


We  are  now  in  a  position  to  derive  the  unique,  stable,  con¬ 
tractive  *  in  the  parametrization  of  all  error  systems  which  max¬ 
imizes  the  entropy  I(E;7;s„). 

Theorem  3.3.  Consider  the  class  of  error  systems  E  which  sat¬ 
isfy  the  condition  ||E[|«,  <7  as  parametrized  in  Lemma  3.1  by 


E  =  7^(R..  +  Qa 


(0  *)>' 


*  £  JIH00,||*||00  <1.  (24' 


Then  the  entropy  I(E;  7;  j„ )  attains  its  maximum  over  this  class 
of  E  with  the  unique  cAoice 


*  =  q  :.(*). 


125) 
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Proof:  Herr  we  adapt  the  proof  in  [2]  to  the  present  setting. 

[  emma  3  i  gives  a!!  r'"r  systems  in  the  form  (24),  where  Raa  + 
is  a.,  pass.  Also, 

detfl  -  [R..  +  Q..]„  (  J  l  ))  =  det( I  -  Q..*) 

which  is  a  unit  in  77H,*,  because  both  Q44  and  $  are  in  TZHoo 
with  HQ..IU  <  1  and  11*11®  -  1  (Lemma  3.1).  Hence  we  may 
apply  Lemma  3.2  to  E  to  obtain 

I(E:  7;  s0 )  =  75 1(*i  li  •’o )  +  7’I([R,ia  +  Qa«]i  i !  1;  •So ) 

-  7’(fle  s0 )/n|def(I  -  Q„(a.)#(s.))|.  (25) 

If  Q4\(s0)  =  0  then 


Proof:  Equation  (31)  follows  immediately  from  equation  (30) 

on  setting  $  =  0.  To  show  (32),  recall  that  R».  +  Q..  is  all 
pass  i.e. , 

(R..  +Q..)’(R..  +  Q..)  =  I.  **  =  }"■  (33) 

The  (1,1)  block  of  this  gives,  Vs  =  j to, 

I  —  [R.a  +  Q..](,  [R«*  +  Q..]i  i 

=  [R«.  +  Q..J),  [R,«  t  Q««]ji  (34) 
so  that,  along  the  imaginary  axis, 

h\det(l-[Ra.  +Q..JMR..  +  Q..]„)| 

=  2ln\det[R..  +Q..;„|  (35) 

=  2M|def(R;i)|  +  2/n|def(Q>J)|,  l36) 


I(E;  7 :  a, )  =  7JI(* ;  1;  s„ )  +  7’I([R..  +  Q„], , ;  1;  a.)  (26) 

which  is  clearly  maximized  by  the  unique  choice  $  =  0  = 
Q),(s0),  and  there  is  nothing  more  to  prove.  So,  henceforth 
in  this  proof  assume  Q, 4(s„)  ^  0. 

Define  the  constant  matrix  H  =  J  by 


where  (36)  follows  from  (35)  on  examination  of  the  structure  of 
Raa  +  Qaa  in  Lemma  3.1 . 

Substituting  (36)  into  the  definition  of  entropy,  we  see  that 


7jI([R..  +  Q..]i  1  i  1;  a0  )  =  —  f  {/njde((R;, 
T  J-oo 


(;^))l 


(  -q  ;.(j.)  (i-q;4(sc}q44(s»))‘/j\ 

+  /n|det(Q4,(>ui))i} 

\  (I  —  Q44(so)Q4’.(s«)),/j  Q44(s0)  ) 

,-371 


(27) 

where  (  •  )1/J  denotes  Hermitian  square  root.  It  is  easy  to  verify 
that  H  is  unitary  and  that 

det(I-H41S(s))  =  def(I-Q4,(s„)S(s)) 
which  is  a  unit  in  . 

Let  us  map  the  unit  ball  in  72H*,  onto  itself  by  the  linear 
fractional  map 


|a0  -  M  J 

Since  R^,  and  Q,,  are  units  in  [8j,  Poisson  s  integral 

theorem  may  be  used  to  evaluate  (37)  as 

7jI([R..  +Q..]ii;1;JIo)  =  l" {Re  sa)  {ln\det{R]^sc))l 

+/n|*f(Q.,(j,))|).  (38) 

The  second  term  in  (31)  is 

-7’ i(Q«\U  );!;*.) 


*  =  ^(H,«),  |f*|U  <  1-  (28) 

Note  that  this  maps  $  =  Q)4(s0)  onto  <£  =  0. 

Lemma  3.2  is  applicable: 

I(*:l:s.)  =I(#;l;s,)  +  I(Q;,(s(,);l;s.) 

-  {Re  So  )ln\det{l  -  Q44(so)*(s,,))|.  (29) 

Lse  this  together  with  (25)  to  relate  the  entropy  of  E  to  the 
entropy  of  *1>: 

I(E;7;  -So )  =  7JI($;  1;  Jo) 

*7’I([R.o  +  Q..)..;1;jo)-7jI(Q;.(j.);1;j.),  (30) 

from  which  it  is  immediate  that  I(E;  7;  s„ )  is  maximized  by  the 
unique  choice  #  =  0.  But  i  —  0  «-»  4>  =  Qj4(s„)  from  above, 
and  the  theorem  is  proved.  AAA 

Denote  maximum  entropy  quantities  by  (  •  )M  E  .  An  expres¬ 
sion  for  the  maximum  entropy  follows  with  ease  from  the  above 
proof. 

Corollary  3.4. 

I(Em  s  ;  7’.  «o ) 

=  73I((R..  +  Q..],,;i;so)  -  yi(Q;.(so);i;s.)  (31) 
=  7 ’(He  s0 ){/n|def(R)1  (s0 ))|  +  fn|*t(Q4,(a0 ))| 

-  ( 1  /2)fn|def(I  -  QI,  (s0  )Q,, (s0 )))}.  (32) 


=  -^/n|det(I-Q,.(so)Q;.(s0))i  /“  U^-]’  i, 

2*  3-oc  1!^  -;w|J 

=  ~^ln\det(l~  Q44(s„  )Q;4(s0  ))|.ir.(/fe  sj 

and  this  with  (38)  gives  (32).  AAA 

III. 2  Entropy  at  infinity 

We  turn  our  attention  in  this  section  to  the  important  spe¬ 
cial  case  which  occurs  when  s,  —*  oo  along  the  real  axis.  The 
entropy  at  infinity  of  G,  (||G||„  <  7)  is  then 

I(G;t;oo )  =  j- J  ln\det(l  -  7*JG,0'ui)G(;w))|  cL;  (39) 

which  is  finite  if  G  is  strictly  proper.  For  our  problem,  where  we 
maximize  the  entropy  of  the  error  system  E,  this  means  that 


the  maximum  entropy  at  infinity  is  finite  if  Ew  E  i  s;rictly 
proper;  this  occurs  when  R  in  the  distance  problem  ,131  is 
strictly  proper,  which  in  turn  occurs  when  P, ,  ( xpl  -  0  the 

standard  configuration  of  figure  1.  This  corresponds  n..  i  rest 
feedthrough  terms  from  the  exogenous  inputs  w  to  1  i'*ri 
outputs  z. 

The  results  for  the  maximum  entropy  problem  ai  it  .fr.it  v  ai  - 
particularly  simple  The  maximum  entropy  solu?<o:  ",  :d 

by  setting  the  arbitrary  stable  contraction  4*  to  .  by 

choosing  the  ‘central  solution’  out  of  the  set  of  poss.blt  K  and 
an  explicit  formula  for  the  maximum  entropy  is  derived  ■  rins 
of  the  state  space  realizations  inherent  in  the  soict:.  ne 

distance  problem  of  Lemma  3.1;  these  state  -pate  re,  ,  .  im  ms 
are  stated  in  the  nex*  lemma 
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Lemma  3.5.  [8]  Consider  the  distance  problem  of  Lemma  3.1 
in  the  cose  R(oo)  =  0.  Suppose  R  has  a  realization 


R  = 

Then  R„a  +  Qaa 
distance  problem 


where 

A:=( 
B  :=  ( 


2* 

"»1  -Pj 

n 

w  /  Lj 

B, 

B, 

"-it  c, 

0 

0 

"'ll  \C, 

0 

0 

(40) 


in  the  parametrization  of  all  solutions  to  the 
||E|U  -  7  ,n  Lemma  3.1  has  a  realization 


R»a  +  Q 


■•-(&$) 


(41) 


D  :  = 


Further, 


A  J 

' 1  /,B, 

7-''»b,  -7-3/,xc;  0 

\ 

0  7- 

J/’Z"YB,  -7  -1/JC[  7-l/1Z- 

•c ;) 

0  \ 

-7-3/,C,X 

■7-3/jb;  Y 

7* '  /3B[Z" 

0 

-7'w,b;  / 

0  I  0' 
0  0  I 
0  0  0 
10  0, 

) 

X  :=  Ric  ( 

'  A-  7-1c;cI\ 

, -B,B(  -  B,B;  -A  J 

(42) 

Y  :=  flic  ( 

A  7-jb,b;\ 

, -C;C,  -c;c,  -A'  ) 

(43) 

Z:=7-’XY-I 

(44) 

A  :=  -A* 

-7-’z--yb3b;  -7-jc;c,x. 

(45) 

Applying  the  results  of  the  previous  section  using  this  realization, 
and  taking  s0  — *  oo  gives  us  the  following  important  theorem. 

Theorem  3.6.  The  entropy  at  infinity,  I(E;7;oo),  is  maximized 
over  the  class  of  error  systems  E  in  the  distance  problem  of 
Lemma  3.1  by  the  unique  choice  =  0. 

If  R(oo)  =  0,  then  the  maximum  entropy  error  system  is 
simply 


E  me  =  7(R««  +  Q«.]n 


/  A 

0 

B, 

Bj 

0 

A 

0 

q-'Z-'YB, 

c, 

0 

0 

0 

Vc, 

— 7*'C,X 

0 

0 

(46) 


and  the  maximum  entropy  is 


I(Eme;  r  oq)  =  7>ace(B|  YB,]  —  7Vace[B'Z~’ YB2)  (47) 
=  7Yace[C, XC’]  —  7Yace[C,XZ_‘C'].  (48) 


Proof:  From  theorem  3.3  the  maximum  entropy  error  system 
is  characterized  by  the  unique  choice  ♦  =  Q',(s0).  Letting  s0  — * 
oo  along  the  real  axis  gives  the  maximum  entropy  at  infinity 
choice  as  ♦  =  Q(,(oo)  =  0,  because  Q,,  is  strictly  proper  from 
Lemma  3.1.  For  details  of  the  limiting  argument  see  [14).  That 
the  maximum  entropy  error  system  has  a  realization  (46)  follows 
easily  by  setting  $  =  0  in  (13)  and  using  Lemma  3.5. 


To  obtain  the  value  of  the  maximum  entropy,  we  take  the 
limit  as  s0  — *  oo  along  the  real  axis  of  the  result  of  corollary  3.4. 
that  is, 

I(E«  7;  oo)  =  lim  (7^0  {/n|«fef(R;i(s<())|  +  /n|det(Q,,(s,  )M 
*0—00 

-(l/2)Hdet(I-Q:,(s0)Q„(So))|}).  (491 

Consider  a  typical  term  from  (49): 

lim  (s0  In  det(I  4-  C(s0I  -  A.)-1  B)). 

*0  — •  OO 

(We  have  dropped  the  modulus  sign  because  sc  is  real  herel. 
Now, 

C(s0I-A)-,B  =  CBs0-1  +  0(s" J ) 
so  by  Lemma  A1  (of  the  appendix)  we  have 

So  In  det{ I  +  C(s„I  -  A)"  ’  B)  =  s„  (Trace[CBs0-  ’ )  +  O(s0'3)) . 
Therefore,  on  taking  limits  as  s„  — *  00, 


lim  (s0  In  det(I  +  C(s0I  —  / B))  =  TracetCB].  (501 

i0—  OO 

Apply  this  to  the  terms  in  (49)  using 

r;,(so)  =  i  +  7'3b;(s,i  + A‘)--YB,  (sn 
and  Q„(so)  =I-7'3B;(soI-A)-,Z--YB3  (521 

from  (41)  to  get 

I(Em e  1 7! 00)  =  7>ace[B[YB;]  -  TracefB'Z-' YB3]  (53) 


as  required,  where  the  third  term  in  (49)  is  zero  in  the  limit 
because  Q,,  is  strictly  proper. 

The  alternative  expression  (48)  follows  in  an  entirely  similar 
fashion;  one  notes  that  I(EM  E ;  7;  00)  =  I(E^  e ;  7;  00)  leading  * 


I(Em  e ;  7;  00)  =  ,*1^(7’  {fn|det(R;a(s„))|  +  /n|det(Q,,(s0 ))( 

-(l/2)fn|det(I-Q.<(so)Q;,(s0))|})  (54) 

which  gives  (48)  in  the  limit.  AAA 

Remark  3.7.  Notice  that  the  entropy  formulae  (47)  and  (48) 
depend  only  on  the  state  space  realization  of  R  and  the  solutions 
X  and  Y  to  the  two  Riccati  equations  (42)  and  (43)  which  are 
inherent  in  the  solution  to  the  distance  problem.  Calculation 
of  the  maximum  entropy  therefore  imposes  negligible  computa¬ 
tional  problems.  Furthermore,  the  maximum  entropy  error  sys¬ 
tem  (46),  being  the  linear  fractional  map  of  $  =  0,  is  simply 
7  times  the  p,  by  m,  (2,1)  block  of  R,,  +  Q..,  which  is  also 
available  from  the  solution  to  the  distance  problem  with  no  extra 
computation. 
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Remark  3.8.  Recall  from  Lemma  2.1  that  { —  I(G;  oc;  oo)}1 11  = 
G,j|  fo'  strictly  proper  G.  Thus  if  we  let  7  — >  00  in  our  maxi¬ 
mum  -'r-tropy  soir  on  we  should  obtain  exactly  the  Hj  optimal 
somtion.  We  show  here  that  this  is  indeed  the  case. 

By  using  the  results  of  Whinner  [  1 3]  it  may  be  shown  that 
the  positive  semidefinite  matrices  -X  and  — Y  are  are  monoton- 
lcaily  decreasing  as  7  increases,  and  that  —  Z  is  monotonically 
increasing  as  7  increases.  Taking -y  — *  oc  we  obtain  ^  =  -I, 


1  with  an  obvious  notation)  which  identifies  the  matrices  —  X-,_  ,, 
and  —  Y,_es  as  the  controllability  and  observability  gramians  of 
R(  —si,  respectively. 

L’sing  this  fact,  a  simple  calculation  shows  that 


I(EW  E  :  oc;  oc)  =  Trace[B[  Y,_  „  B,  ] 

—  Trace  (Bj  (Z,_oo)~*Y,_0CB3)  (57) 

=  -Trace  [[B, ,  B3  j*  [-Y,_„  ][B, ,  B3 ]] 

=  -||R!I!  (58) 


and  that  Q,w  t  =  0. 

It  is  well-known  that  the  Q  £  which  minimizes 


ilElii  = 


R,R+q)| 


(59) 


is  Q  =  QHl  =  0  (the  H3  optimal  solution)  and  in  that  case 
i|EH,i|j  =  ||R||j.  Comparing  this  with  (58)  shows  that  we  have 
I(Ev  e  ;  oc;  oc)  =  — ||EH,  ll] ,  illustrating  the  equivalence  between 
the  maximum  entropy  distance  problem  at  7,  s„  — *  00  and  the 
H,  optimal  distance  problem. 

Note  that  if  7  — >  7,,,,  then  E«  r  -*  E,,i,  the  H*  optimal 
so;ution.  It  is  easy  to  prove  that  |I(EW  e  ;  7;  oo)|  is  monotonically 
decreasing  as  7  increases  from  its  optimal  value;  thus  we  have 
shown  how  7  can  be  used  to  move  from  H<„  optimal  to  H3 
optimal  via  the  maximum  entropy  solutions  for  7,,,  <  7  <  00. 
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Appendix 

The  following  technical  lemma  is  needed. 

Lemma  A1  Let  M  be  a  square  matrix  and  c  >  0.  Then 

(i)  In  det(l— cM)  =  —iT race[M]  +  0(t‘). 

(ii)  —In  det(I  -  dM’M)  >  e3  TraceJM'M]. 

Proof:  Part  (i).  Use  the  Faddeev  formula  (Gantmacher  [7j. 
p88)  to  obtain 

det(l—  eM)  =  1  -  t  7Yace[Mj  +  0(e3) 
and  expand  the  logarithm  of  this  as  a  power  series. 

Part  (ii).  Using  a  well-known  inequality, 

-In  det(l  —  e*M'M)  =  -£/n(l  -e’tr‘(M)) 


=  t3  7race[M-MJ.  AAA 
We  may  now  proceed  with  the  proof  of  Lemma  2  1 
Proof  of  Lemma  2.1:  Part  (i).  By  Lemma  A1  we  can  write 

-i(G;7;  io)  =  ~  j~  |(Tra«:[G*0«)G(jw)]) 

Re  s„  1 3 1 


l-So  -;w| 


-tj:a 


Tract 


G(jtv) 


du>  +  0(7-‘) 
Res „  T 


(s  0  +  1  j 


Therefore, 


r-,:  ,  Re  ] 

G(jw  7— —  .-■■■! 
(s0  +  Jul)  j 


}‘ 


<Li  +  0(7" 3] 


— I(G; 7; s„ )  =  ||G(s)(Re  s„)/(s0  +  s)||]  +  0{-»'3).  («) 

Noting  that  the  0(t~  3 )  terms  are  uonnegativc  we  have 
— I(G;  7;  s„ )  >  ]|G(s)(Rc  s„  )/(s0  -r  j)h; 

whilst 

-I(G.ou;s0)  =  [|G(s)(Re  s0)/(sc  *- 

follows  by  taking  7  — *  00  in  (*). 

Part  (ii)  Firstly  note  that  the  integrands  ni  1*  a*:. arc 
monotonically  increasing  with  j,  and  continuous,  am.  .  .n  ■  bv 
dominated  convergence  both  sides  of  (*)  tend  to  a  limn  r.s  s„ 

00.  Each  side  is  finite  because  ||G|j„  <  7  and  G;o  1  1  0  by 
assumption.  The  result  then  follows  in  a  sirmiai  v.a.  i.<  ..  ,ri,of 

of  Part  (0  7„iA 
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